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INTRODUCTION. 

§ 1. A great many if not all of the problems in mathematics 
may be so formulated that they consist in finding from given 
data the values of certain unknown quantities subject to certain 
conditions. We may distinguish different stages in the solution 
of a problem. The first stage we might say is the proof that the 
quantities sought for really exist, that it is possible to satisfy 
the given conditions or, as the case may be, the proof that it is 
impossible. In the latter case we have done with the problem. 
Take for instance the celebrated question of the squaring of the 
circle. We may in a more generalized form state it thus: Find 
the integral numbers, which are the coefficients of an algebraic 
equation, of which t is one of the roots. Thirty years ago 
Lindemann showed that integral numbers subject to these con- 
ditions do not exist and thus a problem as old almost as 
human history came to an end. Or to give another instance 
take Fermat's problem, for the solution of which the late Mr. 
Wolfskehl, of Darmstadt, has left $25,000 in his will. Find the 
integral numbers x, y, z that satisfy the equation 

x^ + y"* = 2", 

where n is an integral number greater than two. Fermat main- 
tained that it is impossible to satisfy these conditions and he is 
probably right. But as yet it has not been shown. So the 
solution of the problem may or may not end in its first stage. 
In many other cases the first stage of the solution may be so 
easy, that we immediately pass on to the second stage of finding 
methods to calculate the unknown quantities sought for. Or 
even if the first stage of the solution is not so easy, it may be 
expedient to pass on to the second stage. For if we succeed in 
finding methods of calculation that determine the unknown quan- 
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tities, the proof of their existence is included. If on the other 
hand, we do not succeed, then it will be time enough to return 
to the first stage. 

There are not a small number of men who believe the task of 
the mathematician to end here. This, I think, is due to the 
fact that the pure mathematician as a rule is not in the habit of 
pushing his investigation so far as to find something out about the 
real things of this world. He leaves that to the astronomer, to 
the physicist, to the engineer. These men, on the other hand, 
take the greatest interest in the actual numerical values that 
are the outcome of the mathematical methods of calculation. 
They have to carry out the calculation and as soon as they do so, 
the question arises whether they could not get at the same result 
in a shorter way, with less trouble. Suppose the mathematician 
gives them a method of calculation perfectly logical and con- 
clusive but taking 200 years of incessant numerical work to 
complete. They would be justified in thinking that this is not 
much better than no method at all. So there arises a third stage 
of the solution of a mathematical problem in which the object is 
to develop methods for finding the result with as little trouble as 
possible. I maintain that this third stage is just as much a 
chapter of mathematics as the first two stages and it will not do 
to leave it to the astronomer, to the physicist, to the engineer or 
whoever applies mathematical methods, for this reason that 
these men are bent on the results and therefore they will be apt 
to overlook the full generality of the methods they happen to 
hit on, while in the hands of the mathematician the methods 
would be developed from a higher standpoint and their bearing 
on other problems in other scientific inquiries would be more 
likely to receive the proper attention. 

The state of affairs today is such that in a number of cases the 
methods of the engineer or the siu^eyor are not known to the 
astronomer or the physicist, or vice versa, although their prob- 
lems may be mathematically almost identical. It is particularly 
so with graphical methods, that have been invented for definite 
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problems. A more general exposition makes them applicable 
to a vast number of cases that were originally not thought of. 
In this course I shall review the graphical methods from a 
general standpoint, that is, I shall try to formulate and to teach 
them in their most generalized form so as to facilitate their 
application in any problem, with which they are mathematically 
connected.^ The student is advised to do practical exercises. 
Nothing but the repeated application of the methods will give 
him the whole grasp of the subject. For it is not sufficient to 
understand the underlying ideas, it is also necessary to acquire a 
certain facility in applying them. You might as well try to learn 
piano playing only by attending concerts as to learn the 
graphical methods only through lectures. 

^ For the literature of the subject see ''EncyklopSdie der mathematischen 
Wissenschaften/' Art. R. Mehmke, '^Numerisches Rechnen,'' and Art. F 
Willers and C. Runge, "Graphische Integration.'' 
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CHAPTER L 
Graphical Calculation. 

§ 2. Graphical Arithmetic. — ^Any quantity susceptible of mensu- 
ration can be graphically represented by a straight line, the 
length of the line corresponding to the value of the quantity. 
But this is by no means the only possible way. A quantity 
might also be and is sometimes graphically represented by an 
angle or by the length of a curved line or by the area of a square 
or triangle or any other figure or by the anharmonic ratio of four 
points in a straight line or in a variety of other ways. The 
representation by straight lines has some advantages over the 
others, mainly on account of the facility with which the ele- 
mentary mathematical operations can be carried out. 

What is the use of representing quantities on paper? It is a 
convenient way of placing them before our eye, of comparing 
them, of handling them. If pencil and paper were not as cheap 
as they are, or if to draw a line were a long and tedious under- 
taking, or if our eye were not as skillful and expert an assistant, 
graphical methods would lose much of their significance. Or, 
on the other hand, if electric currents or any other measurable 
quantities were as cheaply and conveniently produced in any 
desired degree and added, subtracted, multiplied and divided 
with equal facility, it might be profitable to use them for the 
representation of any other measurable quantities, not so easily 
produced or handled. 

The addition of two positive quantities represented by straight 
lines of given length is eflFected by laying them off in the same 
direction, one behind the other. The direction gives each line a 
beginning and an end. The beginning of the second line has to 
coincide with the end of the first, and the resulting line represent- 
ing the sum of the two runs from the beginning of the first to 
2 1 
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the end of the second. Similarly the subtraction of one positive 
quantity from another is effected by giving the lines opposite direc- 
tions and letting the beginning of the line that is to be subtracted 
coincide with the end of the other. The result of the subtrac- 
tion is represented by the line that runs from the beginning of 
the minuend to the end of the subtrahend. The result is positive 
when this direction coincides with that of the minuend, and nega^ 
tive when it coincides with that of the subtrahend. This leads 
to the representation of positive and negative quantities by lines of 
opposite direction. The subtraction of one positive quantity from 
another may then be looked upon as the addition of a positive and 
a negative quantity. I do not want to dwell on the logical explana- 
tion of this subject, but I want to point out the practical method 
used for adding a large number of positive and negative quantities 
represented by straight lines of opposite direction. Take a 
straight edge, say a piece of paper folded over so as to form a 
straight edge, mark a point on it, and assign one of the two 
directions as the positive one. Lay the edge in succession over 
the different lines and run a pointer along it through an amount 
equal in each case to the length of the line and in the positive 
or negative direction according to the sign of the quantity. The 
pointer is to begin at the point marked. The line running from 
this point to where the pointer stops represents the sum of the 

given quantities. The advan- 
tage of this method is that the 
intermediate positions of the 
pointer need not be marked pro- 
vided only that the pointer keeps 
its position during the move- 
■ ment of the edge from one line 
to the next. As an example take 
the area, Fig. 1. A number of 
rectangular strips ^ cm. wide are substituted for the area so that, 
measured in square centimeters, it is equal to half the sum of 
the lengths of the strips measured in centimeters. The straight 




Fig 1. 
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edge is placed over the strips in succession and the pointer is 
run along them. The edge is supposed to cany a centime- 
ter scale and the pointer is to begin at zero. The final position 
of the. pointer gives half the value of the area in square centi- 
meters. The drawing of the strips may be dispensed with, their 
lengths being estimated, only their width must be shown. K 
the scale should be too short for the whole length, the only thing 
we have to do is to break any of the lengths that range over the 
end of the scale and to count how many times we have gone 
over the whole scale. I have found it convenient to use a little 
pointer of paper fastened on the runner of a slide rule so that it 
can be moved up and down the metrical scale on one side of the 




Fig. 2. 



slide rule. . The area is in this manner determined rapidly and 
with considerable accuracy, very well comparable to the ac- 
curacy of a good planimeter. If the area of any closed curve 
is to be found, the way to proceed is to choose two parallel 
lines that cut off two segments on either side (see Fig. 2), to 
measure the area between them by the method described above 
and to estimate the two segments separately. If the curves of 
the segments may with suflScient accuracy be regarded as arcs 
of parabolas the area would be two thirds the product of length 
and width. If not they would have to be estimated by substitute 
ing a rectangle or a number of rectangles for them. 
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In the same way the addition and subtraction of pure numbers 
may also be carried out. We need only represent the numbers 
by the ratios of the lengths of straight lines to a certain fixed 
line. The ratio of the length of the sum of the lines to the length 
of the fixed lines is equal to the sum of the numbers. The con- 
struction also applies to positive and negative mmibers, if we 
represent them by the ratio of the length of straight lines of 
opposite directions to the length of a fixed line. 

In order to multiply a given quantity c by a given number, 
let the number be given as the ratio of the lengths of two straight 
lines a/6. If the quantity c is also represented by a straight line, 
all we have to do is to find a straight line x whose length is to 
the length of c as a to 6. This can be done in many ways by 





Fig. 3. 
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constructing any triangle with two sides equal to a and b and 
drawing a similar triangle with the side that corresponds to b made 
equal to c. As a rule it is convenient to draw a and b at right 
angles and the similar triangle either with its hypotenuse parallel 
(Fig. 3) or at right angles (Fig. 4) to the hypotenuse of the first 
triangle. Division by a given number is eflFected by the same con- 
struction; for the multiplication by the ratio a/b is equivalent 
to the divisions by the ratio b/a. 

If a, b, c are any given numbers, we can represent them by the 
ratios of three straight lines to a fixed line. Then the ratio of 
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the line constructed in the way shown in Fig. 3 and Fig. 4 to 
the fixed line is equal to the number 

ac 
b' 

Multiplication and division are in this way carried out simul- 
taneously. In order to have multiplication alone, we need only 
make b equal 1 and in order to have division alone, we need only 
make a or c equal 1. 

In order to include the multiplication and division of positive 
and negative numbers we can proceed in the following way. Let 
the lines corresponding to a, x, Fig. 3, be drawn to the right side 
of the vertex to signify positive numbers and to the left side to 
signify negative numbers. Similarly let the lines corresponding 
to b, c be drawn upward to signify positive numbers and down- 
ward to signify negative numbers. Then the drawing of a 
parallel to the hypotenuse of the rectangular triangle a, b through 
the end of the line corresponding to c will always lead to the 

number 

ac 

whatever the signs of a, b, c may be. 

The same definition will not hold for the construction of Fig. 4. 
If the positive direction of the line corresponding to a is to the 
right and the positive direction of the line corresponding to b is 
upwards then the positive directions of x and c ought to be such 
that when the right-angled triangle x, c is turned through an 
angle of 90^ to make the positive direction of x coincident 
with the positive direction of a, the positive direction of c coin- 
cides with the positive direction of 6. If we wish to have the 
positive direction of x upward, the positive direction of c would 
have to be to the left, or if we wish to have the positive direction 
of c to the right, the positive direction of x would have to be 
downward. If this is adhered to, the construction for division 
and multiplication will include the signs. 
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§ 3. Integral Functions. — ^We have shown how to add, subtract, 
multiply, divide given numbers graphically by representing them 
as ratios of the lengths of straight lines to the length of a fixed 
line and finding the result of the operation as the ratio of the 
length of a certain line to the same fixed line. By repeating 
these constructions we are now enabled to find the value of any 
algebraical expression built up by these foiur operations in any 
succession and repetition. Let us see for instance how the values 
of an integral function of x, that is to say, an expression of the form 

Oo + aix + 020:^ + • • • + anPC^^ 

may be found by geometrical construction, where a©, ai • • • an, a: 

are any positive or negative 
numbers. We shall first as- 
sume that all the numbers are 
positive, but there is not the 
least difficulty in extending 
the method to the more gen- 
eral case. 

Now let Oo, di, <h, '"an 
signify straight lines laid off 
on a vertical line that we call 
the j^-axis, one after the other 
as if to find the straight line 
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Oo + c^l 4" 02 + • • • 4" On. 



The lengths of these lines measured in a conveniently chosen 
unit of length are equal to the numbers designated by the same 
letters. In Fig. 5 do runs from the point to point Ci, ai from 
Ci to ft, • • • On from C« to C„+i. 

Let a: be the ratio of the lines Ox and 01, Fig. 5, drawn hori- 
zontally from to the right. The length 01 is chosen of con- 
venient size independent of the unit of length that measures the 
lines Oo, oi, • • • o». The length Ox is then defined by the value 
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of the ratio x. Through x and 1 draw lines parallel to the 2^-axis. 
Through Cn+i draw a line parallel to Oxj that intersects those two 
parallels in Pn and Bn> Draw the line BnCn that intersects the 
parallel through x in Pn-i. Then the height of P»-i above (7« 
will be equal to On^. For if we draw a line through Pf»-i parallel 
to Ox intersecting the y-axis in Z)„, the triangle CJDnPn^i will be 
similar to CnCn^iBn and their ratio is equal to x, therefore 
CnDn = anX. Couscquently the height of P«_i above Cnr-i is 
equal to C»-iZ)» = anX + a^^i. Now let us repeat the same 
operation in letting the point Dn take the part of Cn+i. Through 
Dn draw a line parallel to Ox, that intersects the parallels through 
X and 1 in P»-i and Bnn-i> Draw the line Bn^iCn-i that intersects 
the parallel through x in Pn^2. 
Then the height of P»-2 above 
Cn^i will be equal to 

Cn^iDn'X = (anX + a„^i)x, 
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and the height above C»-2 will be 
equal to 

Continue in the same way. Draw 
Pnr-%Bnr-2 parallel to Ox, draw 
Bnr^2Cnr^ and find the point P«^. 
Then the height of P»-8 above Cn-^ will be 

and the height of P»-8 above C«^ 

Finally a point Po is found (see Fig. 6 f orn = 4) by the inter- 
section of BiCi with the parallel to the y-axb through x, whose 
height above is equal to 

anX"^ + an-i**"* + • • • + aia + oo 
Let us designate the line xP^ by y, so that 
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in the sense that y is a vertical line of the same direction and 

length as the sum of the vertical lines a^ix;'^, an^io;"^^, • • • aix, a©. 

The same construction holds good for values of x greater than 

1 or negative. The only difference is that the point x is beyond 

the interval 01 to the right of 1 or to the left of 0. The negative 

sign of 

anXf anX + an-i, anX^ + o»-i«, etc., 

will signify that the direction of the lines is downward. Nor are 
any alterations necessary in order to include the case that several 
or all of the lines a©, ai, • • • an are directed downward and corre- 
spond to negative numbers. They are laid off on the j^-axis in 
the same way as if to find the sum 



cto + ^1 "i" ^2 + • • • + a 



n> 



0«+i lying above or below C. according to a. being directed 
upward or downward. The construction can be repeated for a 
number of values of x. The points Po will then represent the 
ciurve, whose equation is 

y = ao + aia; + h a^i^'^, 

x and y measuring abscissa and ordinates in independent units 
of length. 

In order to draw the curve for large values of a; a modification 
must be introduced. It will not do to choose 01 small in order 
to keep X on your drawing board; for then the lines 5.0. will 
become too short and thus their direction will be badly defined. 
The way to proceed is to change the variable. Write for instance 
X ~ :r/10, so that X is ten times as small as x and write 

^.= a..lO*. 
Then as 

X XT X 

y = ao + ai-10jQ+a2lO^^H hanlO^j^ 
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we find 

Lay off the lines Ao, Ai, - • • An in a convenient scale and let 
X play the part that x played before. The curve differs in scale 
from the first curve and the reduction of scale may be different 
for abscissas and ordinates but may if we choose be made the 
same so that it is geometrically similar to the first curve reduced 
to one tenth. It is evident that any other reduction can be 
effected in the same manner. By increasing the ratio x/X we 
enhance the value of An in comparison to the coefficients of lesser 
index, so that for the figure of the ciu^e drawn in a very small 
scale all the terms will be insignificant except AnX^. In this 
case the points Ci, C2, • • • , Cn will very nearly coincide with 
and only CnH will stand out. 

It is interesting to observe that the best way of calculating an 
integral function 

(lo + Oia: + (hix^ + • • • + dnX^ 

for any value of x proceeds on exactly the same lines as the 
geometrical construction. The coefficient an is first multiplied 
with X and a„r-i is added Call the result an^i. This is again 
multiplied by x and a„^2 is added. Call this result a»_2'. Con- 
tinuing in this way we finally obtain a value of ao', which is equal 
to the value of the integral function for the value of x considered. 
Using a slide rule all the multiplications with x can be effected 
with a single setting of the instrument. The coefficients a^ and 
the values a/ are best written in rows in this way 

an anr~i a»_2 * * * di ao 
anX ttn^ix • • • a^'x aix 



On-X anr-2 ' " tti tto 



The accuracy of the slide rule is very nearly the same as the 
accuracy of a good drawing. But the rapidity is very much 
greater. When therefore only a few values of the integral func- 
tion are required, the geometrical construction will not repay 
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the trouble. It is diflFerent, however, when the object is to make 
a drawing of the curve. The values supplied by calculation 
would have to be plotted, while the geometrical construction 
furnishes the. points of the curve right away and in this manner 
gains on the numerical method. 

There is another geometrical method, which in some cases 
may be just as good. Let us propose to find the value of an 
integral function of the fourth degree. 

y =^ ao + aio; + ani? + a^pf + a^ 

and let all coefficients in the first instance be positive. 

The coefficients Oo, ai, 029 os, a^ are supposed to be represented 
by straight lines, while x will be the ratio of two lines. The lines 

ao, ai, 02, os, a^ are laid off in a 
broken line ao to the right from 
Co to Ci, ai upward from Ci to 
C2, 02 to the left from C2 to C3, Os 
downward from Cz to C4, 04 again 
to the right from (74 to Cs (Fig. 7). 
Through Cs draw a line (7^4 to 
a point A on CzC^ or its prolonga- 
tion and let x be equal to the 
ratio C4A : C^Ph taken positive 
when ^4^ has the same direc- 
tion as CzC^. Then we have 




and 



C4A = (kPCf 
CzA = 042: + Os. 



CiA and dzA are positive or negative according to their direction, 
being the same as the direction of CzCa or opposite to it. Through 
A draw the line AB forming a right angle with C^ to a point B 
on CiCs or its prolongation. Then we have 



and 



CzB = CzA'X = (04a: + (h) X 



GRAPHICAL CALCDIATION. 11 

CtB = o^a:* + ojl + Oi. 

CtB and dB are positive or negative according to thrar direction 
being the same as the direction of dCt or opposite to it, Snii> 
larly we get 

CiD = 0*3^ + Ota? + oja; + Ob 
and finally 

CoE = 0421* + OiS? + oja? + aix + oq. 

CoE U positive, when E is on the right side of Co and negative 

when on the left side. When the point A moves along the line 

CiCi, the point E will move 

along the line CoCi and its 

position will determine the 

values of the integral function. 

To find the portion of E for 

any position of A, we might 

use transparent squared paper, 

that we pin onto the drawing 

at Cfe so that it can freely be 

turned round Cg. Followine 
..... t Fio. 8. 

the Imes of the squared paper 

along ChABDE after turning it through a small angle furnishes 

the position of E for a new portion of A (Fig, 8). 

To include the case of negative coefficients we draw the corre- 
sponding line in the opposite direction. If for instance Oj is 
negative CtCi would have to lie above C3; but CtA would have 
to be counted in the same way as before, positive in a downward, 
negative in an upward direction. 

The extension of the method to integral functions of any degree 
is obvious and need not be in^sted on. It may be applied with 
advantage to find the real roots of an equation of any degree. 
For this purpose the broken line CsABDE would have to be 
drawn in such a way that E coincides with Co. lu the case of 
Fig. 7, for instance, it is easily seen that no real root exists. 
Fig. 9 shows the application to the quadratic equation. A circle 
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is drawn over CoCs as diameter. Its intersections with CiCt 
furnish the points A and A' that correspond to the two roots. 
Both roots are negative in this case. 

The first method of constructing 
the values of an integral function can 
be extended to the case where the 
fimction is given as the sum of a 
number of polynomials of the form 



/ 



A1 



^-#= 



^^^ 



Oo 



y = oo + ai(x - p) + (h{x—p)(x — g) 
+ az(x — p)(x — q)(x — r) + • * •• 



Ql 
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Let us again suppose oo, at, 02, • • • 
to represent straight lines laid off as 
before on the y-axis upwards or down- 
wards as if to find their sum. a:,p,g,r ... are meant, to be num- 
bers represented by the ratio of certain segments on the axis of 
abscissas. Let us consider the case of four terms, the highest poly- 
nomial being of the third degree. The fixed distance between the 
points marked p and p + 1^ 9 and q+ 1, r and r + 1 on the 
ads of abscissas. Fig. 10 is chosen arbitrarily and the position 



a. 






^:C. ^- 




a 



X p^\ q¥\ ^*L 
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of the points marked p, q, r, x is made such that the ratio of 
Op, Oq, Or, Ox to that fixed distance is equal to the numbers 
p, q, r, X. For negative values the points are taken on the left 
ofO. 



GBAPmCAL CALCULATION. 13 

Draw parallels to the y-axis through j>, q, r, x, p+ 1, q+ 1, 
r + 1. On the parallel through r + 1 find the point Qo of the 
same ordinate as d and on the parallel through r find the point 
Ao of the same ordinate as Cz. Join Aq and Qo by a straight 
line and find its intersection Pi or that of its prolongation with 
the parallel through x. The height of Pi above Cz or Aq is 
equal to az(x — r) and the height above C2 is equal to 03(0; — r) 
+ 02. On the parallel through q+ I find a point Qi of the same 
ordinate as Pi and on the parallel through q a point Ai of the 
same ordinate as C2. Join Ai and Qi by a straight line and find 
its intersection P2 or that of its prolongation with the parallel 
through X. The height of P2 above C2 or Ai is equal to 

[az{x — r) + a2](x— q), 

and the height above Ci is equal to 

az(x — r)(x — q) + (hix — q) + (h. 

Finally find a point Q2 on the parallel through p + 1 of the 
same ordinate as P2 and a point A2 on the parallel through p of 
the same ordinate as Ci. Join A2 and Q2 by a straight line and 
find its intersection Pz or that of its prolongation with the par- 
allel through X. The height of Pz above Ci or A2 will then be 

equal to 

[az(x — r){x — q) + (h{x — g) + ai]{x — p) 

and the ordinate of Pz will be equal to the given integral function 

y = a^{x - r){x - q)(x - p) + (h(x - q)(x — p) 

+ ai(x — p) + ao. 

For large numbers p, q, r, x we use a similar device as before by 
introducing new numbers P, Q, R, X equal to one tenths or one 
hundredth or any other fraction of pgraj. For instance 

P=p/10, Q = g/10, fl=r/10 Z=r/10. 

We then write 

Ao =00, ^1 = lOoi, A2 = IOO02, Az = lOOOos, 
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and obtain 

y ^ Ao + Ai(X - P) + A2(X - P)(X ^ Q) 

+ ^(Z"-P)(Z-.Q)(Z-B). 

The scale for the lines Ao, Ai, A2, Az and y must then be reduced 
conveniently and the values are constructed in the same way as 
before. 

Now let us consider the inverse problem. The values of the 
integral function are given for 

x== P, q, r, s; 

find the lines oq, a\, di, a^, so that the value of the integral function 
may be found for any other value of x in the way shown above. 

Let us designate the given values of the integral function for 
^= Py9>^9^ by yp, yq, yr, y. and the points on the parallels through 
p, q, f, 8 with these ordinates by P, Q, R, S (see Fig. 12). 

For a; = p the integral function 

y- Oq + ai{x — p) + 02(0: — p)(x — ?) + (h(x— p){x—q){x'-r) 

reduces to Oo- Therefore we have yp = oq. The point C\ is 
found by drawing a parallel to the axis of abscissas through P 

and taking its intersection with 
the axis of ordinates. 

In order to find C2 draw a 
straight line through P and Q 
and find its intersection A with 
the parallel through p + 1 (Fig. 
11). A parallel to the axis of 
abscissas through A intersects 
the axis of ordinates in C2. For 
the differences yq—yp and ya-^yp 
(writing ya for the ordinate of 
A) are proportional to the differences of the abscissas and con- 
sequently in the ratio (? — p) : 1. Therefore 

yq — yp 

y^^vp^ ~ — ~ = ^1- 
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In the same way as the point Q on the pariJIel through q we 
might join any point Z on a paralld throi^h x with the point P, 
find the intersection with the parallel through p +, 1 and draw a 
parallel to the axis of abscissas. The point of intersection of 



this parallel with the vertical through x let us call X' and its 
ordinate y'. Then we have 

y'~yp = T^r^= «> + «s(* - «) + "^i^ - ?)(« - O. 

X y 
Let us carry out this construction not only for a: = 5 but also 
for a: = r and x = a. This leads us to three points Q', R', S' 
on the verticals through q, r, s, whose ordinates are the values 
of the integral functions 

y' = (flo + fli) + OaCa: — ?) + <h{x — q)(_x — r). 
In this way we have reduced our problem. Instead of having 
to find an integral function of the third degree from four given 
points P, Q, R, S, we have now only to find an integral function 
of the second degree from three given points Q', R', S'. A second 
reduction is effected in exactly the same manner. Q' is joined 
with R' and S' by striught lines and through their intersection 
with the vertical through 9 + 1 parallels to the axis of abscissas 
are drawn that intersect the verticals through r and s in the 
points R" and S" respectively. The ordinates of these points 
are the values of the integral function y" defined by 
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f or a; = f and a; = *, or 

y" = ao + oi + 02 + c^(x — r). 

The horizontal through R" intersects the axis of ordinates in the 
point Cz. Finally we find Ca by drawing a parallel to the axis 
of abscissas through the intersection of JR"iS" or its prolongation 
with the vertical through r + 1. 

Having found the points CiC^CzC^ we can now for any value 
of X construct the ordinate 

y = ao + ai{x — p) + 02(0; — p){x — q) 

+ 03(0; — p){x — q)(x — r), 

and thus draw the parabola of the third degree passing through 
the four points P, Q, iJ, 8. 

The construction may be somewhat simplified first by making 
p + 1 = g. Our data are the points P, Q, R, 8, and we are 
perfectly at liberty to make the vertical through p + 1 coincide 
with the vertical through Q, In this case the point Q' will 
coincide with Q. The parabola of the second degree through the 
points Q'R'8' is again independent of the distance between the 
verticals through q and q+ 1 and at the same time independent 
of the point P. Therefore we are perfectly at liberty, for the 
construction of any point of this parabola, to make the vertical 
through q+ 1 coincide with the vertical through jR even if the 
distance of the verticals through P and Q is different from that 
of the verticals through Q and R. R" will in this case coincide 
with R'. The procedure is shown in Fig. 12. Starting from 
the points P, Q, R, 8 the first step is to find R\ 8' by connecting 
JB and 8 with P and drawing horizontals through the inter- 
sections Ar, At with the vertical through q. The next step is to 
find S" by connecting Q (identical with Q') with 8' and drawing 
a horizontal through the intersection with the vertical through r. 
Now the straight line /J"fii" can be drawn (R" being identical 



GRAPHICAL CALCULATION, 



17 



with R'). On the vertical through any point x take the inter- 
section with U"i8" and pass horizontally to the point Ax' on the 
vertical through r. Draw the line Q'Ax' and find its intersection 
with the vertical through x. This point is on the parabola 
through Q'R'S\ Pass horizontally to the point A^ on the 
vertical through q and draw the line AxP. Its intersection with 
the vertical through x is a point on the parabola of the third 
degree through P, Q, R, S. 

The method is evidently applicable to any number of given 
points, the degree of the parabola being one unit less than the 
number of points. 

The methods for the construction of the values of an integral 
function may be applied to find the value of any rational f imction 

y = R{x). 

For a rational function can always be reduced to the form of a 
quotient of two integral functions 

R(^) = gi(x)lg2(x). 

Now after having constructed curves whose ordinates give the 
values of gi(x) and g2{x) for any abscissa x (Fig. 13), R(x) is found 
in the following manner. 
Through a point P on the 
axis of abscissas draw a 
parallel to the axis of or- 
dinates. Let Oi and O2 
be the points whose ordi- 
nates are equal to gi(x) 
and gi{x). Pass horizon- 
tally from Gi to G\ on the 
vertical through P and 
from 6r2 to G% on the axis of ordinates. Draw a line through 
P and G2 and produce it as far as ^ where it intersects the 
horizontal through 6^. Then R{x) is equal to the ratio GiA 
to PO. Gi'A may then be set off as ordinate on the vertical 
3 
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through X and defines the point M whose ordinate is equal to 
R{x) in lengthy when OP is chosen as the unit of length. 

§4. Linear Functions of Any Number of Variables. — ^Let us 
consider a linear function of a number of variables a;i, o:^ • • • Xn, 

Oo + aiXi + 0^X2 + . . . + a^Xnf 

where Oo, auOii, ••• On are given numbers positive or negative. 
The question is how the value of this linear function may be 
conveniently constructed for various systems xu X2, ••• Xn> 
Suppose ao, ai, • • • an to represent horizontal lines directed to 
the right or left according to the sign of the corresponding number 
and to be laid off on an horizontal axis in succession as if to find 
the sum 

Oo + ai + 02 + • • • + On, 

ao begins at and runs to Ci, a^ begins at Ci and runs to C2 and 
so on (Fig. 14). The numbers xi, X2, • • • Xn let us represent 
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by ratios of lengths. We draw a vertical line through and 
choose a point P on the horizontal axis. Then let xi be equal 
to the ratio 01/PO, X2 = 02/PO, etc. If P is chosen on the left 
of 0, we take the point 1 above for a positive value of xi and 
below for a negative one and the same for the other points. 
Mark a point above in the same distance from as P. Join 
the point P with the points 0, 1, 2, 3, 4, • • • and draw a broken 
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line OA0A1A2AZA1 in such a manner that Aq is on the vertical 
through Ci and OAq is parallel to PO, Ai on the vertical through 
C2 and AqAi parallel to PI, A2 on the vertical through Cb and A1A2 
parallel to P2 and so on. Then the ordinate yo of Ao will have 
the same length as ao and will be directed upward when the 
direction of oo is to the right, and downward when the direction 
of ao is to the left. The difference yi — yo of the ordinates of Ai 
and Ao is equal in length to aiXu as yi — yo and at have the same 
ratio as 01 and PO. Ai will be above or below Aq according to 
the line aiXi being directed to the right or to the left and it is 
understood that aiZi has the same direction as Oi for positive 
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values of xi and a direction opposite to ai for negative values 
of xi. Thus the ordinate yi has the same length as the line 
ao + aiXi and its direction is upward or downward according to 
the direction of the line ao + aiXi being to the right or to the left. 
In the same way it is shown that the ordinate y^ of the point A2 
is equal in length to 

ao + CLiXi + 02252, 

and ys to 

Oo + (llXi + 02252 + 0»358 

and so on, the direction upward or downward corresponding 
to the positive or negative value of the linear function. 
If the values of xi, 0D2, • • • Xn satisfy the equation 

Oo + aiXi + 020:^ + • • ' + dnJXn = 

the ordinate yn must vanish, that is to say, the point An must 



20 GRAPHICAL METHODS. ^ 

coincide with Cn+i> the end of the line an? And vice versa if An 
and Cn+i coincide the equation is satined. Consequently if we 
know all the values but one of the niunbers xi, Xi, • • • Xn the 
unknown value can be found graphically. For suppose xz to be 
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the unknown value we can, beginning from 0, find the broken 
line as far as A% and beginning from the other end An we 
can find it as far as Az (Fig. 15). A parallel to ^42^8 through P 
furnishes the point 3 on the axis of ordinates. Ji xuXi, • • • Xn^i 
are known and only Xn not, we can draw the broken line as far 
as An-\ and as An has to coincide with Cn+i> we can draw a parallel 
to Anr-iAn through P and find the point n on the axis of ordinates 
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that determines the value Xn by the ratio On/PO or On/Oo. In 
Figs. 15 and 16 all the coefficients oq, ai, ">, are positive. A 
negative coefficient as is shown in Fig. 17. The only difference 
is that Cfi lies to the left of Cs and consequently the broken line 
passes from Aa back to A^. 

If we keep the points 0, 1, 2, • • •, in their positions but change 
the position of P to P' (Fig. 18) and repeat the construction of 
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the broken line, we obtain OA^A\Ai • • • instead of OA^\At 
The ordinate y.' of the point A^ is evidently 



00 



01 



y« =• <»o D//1 + 01-577^ + 



fo 



and therefore 



PO 



• • 






y« "~ pfQ y«* 



That is to say, by changing the position of P without changing 
the position of the points 0, 1, 2, • • • we can change the scale of 
the ordinates of the broken line. They change inversely pro- 
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portional to PO. It may be convenient to make use of this 
device in order to make the ordinates a convenient size inde- 
pendent of the scale that we have chosen for the points 0, 1, 2, • • • 
that determine the values 



X\ = 



01 



02 



00' ^""00* 



A linear equation with only one unknown quantity 

Oo "f" tti^^i = 

is solved by drawing a parallel to A^\ through P. Let a second 
equation be given with two unknown quantities 

fco + hxi + hiXi = 0. 

The lines bo, &i, 5^ are laid off as before. Knowing Xi as the 
solution of the first equation we can construct the broken line 
OBqBi corresponding to the second equation and as B2 must 
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coincide with the end of bz, we can draw a parallel to B1B2 
through P and find Xz. In a similar manner we can find xz 
from a third equation 

Co+ C1X1+ O2X2 + 02^ = 0, 

and so we can find any number of unknown quantities, if 
each equation contains one unknown quantity more than those 
before. 

In the general case when n unknown quantities are to be 
determined from n linear equations each equation will contain 
all the imknown quantities, and therefore we cannot find them 
one after the other as in the case just treated. But it can be 
shown that by means of very simple constructions the general case 
is reduced to a set of equations, such as has just been treated. 

Let us begin with two 
equations and two un- 
known quantities. 

do "4" cti3?i ~l~ (W2 = 0, 

The lines Oo, a^, oi are laid oflp on a horizontal line OAqAiAz and 
the lines bo, bi, 62 on another horizontal line O'BqBiBz (Fig. 19). 
Now let us join and 0', Ao and Bo, A\ and Si, A2 and B2 by 
straight lines and let us draw a third horizontal line intersecting 
them in the points C'CoCiCi. These points correspond to a 
certain linear function 

Co + ClXi + C2XZ9 

and it can be shown that it vanishes when x\ and x^ are the same 
values for which the first two linear functions vanish. Let the 
distance of the first two horizontal lines be I and the distance of 
the third from the first and second h and h. Then it can readily 
be seen that 

Co = Oo + "7 (60 "" flto) = T do + "T 6o« 
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For a parallel to 00' through Aq defines with the line AoBq on 
the third and second horizontal line segments equal to Co — oq 
and 60 — do and as these segments have the ratio h/l, it follows 
that 

Co = oo + y (5o — Oo) = 7 ^ + T ^o* 

By drawing a parallel to AqBo through Ai and to AiBi through 
A2 or through Bt (which comes to the same thing), we convince 
ourselves in the same way that 

Ci = ai + y (5i — ai) = y ai + y 6i 
and 

Multiplying the equation 

by k/l and the equation 

60 + 6iiPi + h^ = 
by h/l and adding the two products, we obtain 

Co + CiXi + OiX2 = 0. 

The third horizontal need not lie between the first two. If it 
lies below the second we have merely to give k a negative value 
and if it lies above the first we have to give h a negative value 
and the same formulae for Cq, ci, 02 hold good. Consequently the 
conclusion remains valid, that from the first two equations the 
third follows. 

Now as we are perfectly at liberty to draw the third horizontal 
line where we please, we can let it run through the intersection 
of the straight lines AiBi and A2B2. In this case the points Ci 
and C2 must coincide and consequently 02 must vanish. If ci 
does not vanish we can by what has been shown above find xi 
and with 0:1 we can find X2 from either of the two first horizontal 
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lines. In case ci also vanishes^ that is to say^ in case the three 

straight lines A2B29 AiBi, AqBq all pass through the same pointy 

while 00' does not pass through it, the two given equations 

cannot simultaneously be satisfied. For if they were, it would 

follow that 

Co + Ci^i + O2OD2 = 0, 

and as Ci and Ci are zero cq would have to be zero, which it is not 
as 00' is supposed not to pass through the intersection of .^2-82, 
AiBi and AqEq. If on the other hand all four lines ^2^2^ AiBi, 
AqBo, 00' pass through the same point, cq, ci and c^ will all three 
vanish. In this case the two given equations do not contradict 
one another, but &o&i^ will be proportional to aoaia2. The 
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Fig. 20. 

second equation will therefore contain the same relation between 
xi and QD2 as the first, so that there is only one condition for Zi 
and 0^ to be satisfied. We may then assign any arbitrary value 
to one of them and determine the value of the other to satisfy the 
equation. 

In the case of two linear equations of any number of quantities 
Xi, a^, • • • OTn we can by the same graphical method eliminate one 
of the quantities. In Fig. 20 this is shown for two linear equa- 
tions with six unknown quantities. The two horizontal lines 
OAoAiA2AzAiAsAQ and O'BJBABzBiBJBt represent two linear 
equations. Through the intersection of AzBz and A^Pa^ a third 
horizontal line is drawn intersecting the lines 00', AqBq, AiBi, 
• • • ^6^6 in 0"CoCi '" C%. As Cz and d coincide, the line c* 
vanishes and x^ is eliminated, so that the equation assumes the 
form 

Co + CiXl + OiOD2+CzPCz + c^ixb + c^pc^ = 0. 
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Suppose now that a set of six equations with six unknown quan- 
tities is represented geometrically on six horizontal lines. We shall 
keep one of these; but instead of the other five we construct five 
new ones from which one of the unknown quantities has been 
eliminated by means of the first equation. Now it may happen 
that at the same time another unknown quantity is eliminated^ 
then this quantity remains arbitrary. Of the five new equations 
we again keep one that contains another unknown quantity and 
replace the four others again by four new ones from which this 
unknown quantity has been eliminated. Groing on in this 
manner the general rule will be that with each step only one 
quantity is eliminated, so that at last one equation with one un- 
known quantity remains. Instead of the given six equations 
with six unknown quantities each, we now have one with six, 
one with five and so on down to one with one. The geometrical 
construction shows that this system is equivalent to the given 
system, for we can just as well pass back again to the given 
system. We have seen above how the unknown quantities 
may now be found geometrically. It may however happen in 
special cases that with the elimination of one unknown quantity 
another is eliminated at the same time. To this we may then 
assign an arbitrary value without interfering with the possibility 
of the solution. Finally all unknown quantities may be elimi- 
nated from an equation. If in this case there remains a term 
different from zero it shows that it is impossible to satisfy the 
given equations simultaneously. If no term remains, the two 
equations from which the elimination takes its origin contain the 
same relation between the unknown quantities and one of them 
may be ignored. 

§ 5. The Graphical Handling of Complex Numbers. — ^A complex 
number 

z = z + yi 

is represented graphically by a point Z whose rectangular coordi- 
nates correspond to the numbers x and y. The units by which 
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the coordinates are measured^ we assume to be of equal lengtlu 
We might also say that a complex number is nothing but an 
algebraical form of writing down the coordinates of a point in a 
plane. And the calculations with complex numbers stand for 
certain geometrical operations with the points which correspond 
to them. 
By the "sum" of two complex numbers 

zi = xi + yii and Zi = X2 + y^ 
we understand the complex number 



where 

and we write 



xz- xi + Xi and yz^ yi + jfi, 

258 = Zl + 2f2. 



Graphically we obtain the point Zz representing zz from the 
points Zi and Z2 representing zi and 22 by drawing a parallel 

to OZ2 through Zi and making ZiP 
(Fig. 21) equal to OZ2 in length 
and direction or by drawing a paral- 
lei through Z2 and making Z2P 
equal to OZi in length and direc- 
tion. The coordinates of P are 
evidently equal to X1 + QD2 and 

yi + 2^2. 
Two complex numbers z and zf 

are called opposite, when their sum 

is zero. 




Fig. 21. 



2 + a' = or a: = — a;' and y =^ -- y^ or z = — z'. 

The corresponding points Z and Z' are at the same distance from 
the origin but in opposite directions. 

The difference of two complex numbers is that complex 
number, which added to the subtrahend gives the minuend 

22 + (zi ~ 252) = Zi. 
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Therefore 

zi — 252 = (id — a^) + (yi — y2)i. 

This may also be written 

zi + 22' where 252'= — 2^2 = — a^i — ^t* 

That is to say, the subtraction of the complex number 22 from 211 
may be effected by adding the opposite number — 212. For the 
geometrical construction of the point Z corresponding to 2:1 — 2!^ 
we have to draw a parallel to OZ2 through Zi and from Zi in 
the direction from Z2 to we have to lay off the distance Z2O. 
Or we may also draw from a line equal in direction and in length 
to Z2Z1. This will also lead to the point Z representing the 
difference Zi — 212. 

The rules for multiplication and division of complex numbers 
are best stated by introducing polar coordinates. Let r be the 
positive number measuring the distance OZ in the same unit 
of length in which x and y measure the abscissa and ordinate, so 

that 

r= ^^x^ + f 

and let fp be the angle between OZ and the axis of x, counted in 
the direction from the positive axis of x toward the positive 
axis of y through the entire 
circumference (Fig. 22). Then 

we have 
a; = r cos ^, y = rsui<p 

and 

z= x + yi= r(cos^+ sin <pi). 

Let us call r the modulus 
and <p the angle of 2;. The an- 
gle may be increased or di- 
minished by any multiple of four right angles without altering 
z, but any alteration of r necessarily implies an alteration of z. 
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According to Moivre's theorem, we can write 

z = re'^. 

By the product of two complex numbers 

2i = fi6*" and 2^ = fi^^ 

we understand that complex number Zg whose modulus rz is 
equal to the product of the moduli fi and r2 and whose angle <ps 
is the sum of the angles ^ and ^ or differs from the smn only by 
a multiple of four right angles 

The definition of division follows from that of multiplica^ 
tion. The quotient zi divided by z^ is that complex number, 
which multiplied by z^ gives zi. Therefore the product of its 
modulus with the modulus of z^ must be equal to the modulus of 
zi and the sum of its angle with the angle of z^ must be equal to 
the angle of Zi, Or we may also say the modulus of the quotient 
21/22 is equal to the quotient of the moduli ri/r2 and its angle is 
equal to the difference of the angles ^ — ^. An addition or 
subtraction of a multiple of four right angles we shall leave out 
of consideration as it does not affect the complex number nor 
thepointrepresentingit. 

The geometrical construction corresponding to the multi- 
plication and division of complex numbers is best described by 
considering two quotients each of two complex nimibers that 
give the same result. Let us write 

21/22 = Zz/Zi. 

The geometrical meaning of this is that 

Ti/r2 = TzIUy 
and 

^ — ^ = ^ — ^4. 

That is to say, the triangles ZiOZ% and ZzOZ^ are geometrically 
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similar (Fig. 23). When three of the points Zu Zt, Zt, Zi are 
given the fourth can evidently be found. For instance let 
Zi, Z2, Z4 be given. Draw a parallel to Z1Z2 intersecting OZt 
at a distance u from 0. This point together with the inter- 
section on OZi and with will form the three comers of a tri- 
angle congruent to the triangle Z^^zO. It will be brought into 
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the position of ZiZzO by being turned round so as to bring the 
direction of the side in OZ2 into the position of OZ4. Thus the 
direction of OZ3 and its length may be found. 

This construction contains multiplication as well as division as 
special cases. Let Z4 coincide with the point a; = 1, y = 0, so 
that Z4 = 1 (Fig. 24), then we have 

Zi/Sfz = 28 or 2i = 2f22i8. 

From any two of the points Zi, Z2, Zz a simple construction gives 
us the third. 

The geometrical representation of complex numbers may be used 
to advantage to show the properties of harmonic oscillations. 

Let a point P move on the axis of x^ so that its abscissa at the 
time t is given by the formula 

a: = r cos {rd + a), 

n, r and a being constants. We call r the amplitude and n^ + a 
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the phase of the motion. The point P moves backwards and 
forwards between the limits x= r and a; = — r . The time 
T = 2T/n is called the period of the oscillation, it is the time in 
which one complete oscillation backwards and forwards is per- 
formed. 
Now instead of x let us consider the complex number 

2 = f cos(n^ + a) + r sin(nt + a)i 
or 

of which X is the abscissa and let us follow the movement of 
the point Z. For < == we have 

Designating this value by Zo, we can write 
The geometrical meaning of the product 

is that the line OZq is turned round through the angle rU. For 
the modulus of e"'* being equal to 1 the modulus of zo is not 

changed by the multiplication. The 

movement of the point Z therefore 

•n. consists in a uniform revolution of 

^t«o ^^ round 0. At the moment t= 

,01-'"^. the position is OZq and after the 

^'^ j > time T = 2x/n the same position is 

\ / occupied again. The revolution goes 

\ / on in the direction from the positive 

— -^ ^"^*^ axis of aj to the positive axis of y 

(Fig. 25). 
j^Q 2g^ The movement of Z is evidently 

simpler than the movement of the 
projection P of Z on the axis of x. 
Let us consider a motion composed of the sum of two harmonic 



1 — 



t-T^' 



1 
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motions of the same period but of different amplitudes and 

phases 

a; = fi cos (rrf + oi) + f2 cos (frf + as), 

and let us again substitute the motion of the point Z correspond- 
ing to the complex number 

For < = the first term is 
and the second term 

Introducing 21 and s^ into the expression for z we have 

z = 2i6**** + afec*** == (zi + sfe)e**' = a^*** 
where 

2Sb = 2l + 2S8. 

This shows at once that the movement of Z is a uniform circular 
movement consisting in a uniform revolution of OZ round 0. 
The position at the moment < = is OZi corresponding to the 
complex number 

2b = 2i + %• 

The projection of Z on the axis of x has the abscissa 

a; = fs cos (n< + as) 

where fs and as designate modulus and angle of 7^ Thus the 
sum of two harmonic motions of the same period is shown also 
to form a harmonic motion. 

The same holds for a sum of any number of harmonic motions 
of the same period. For the complex number 

where ri, r%, • • • fx; ai, a^, •*• a^ and n are constants may be 

written 

z = zie^ + ZiB*** + • • • + Zifi^ 
or 
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where 

a© =* zi + 2i + • • • + Zx' 

The movement of Z therefore^ excepting the ease zo = 0, consists 
in a uniform revolution of OZ round 0, OZ always keeping the 
same length equal to the modulus of zq. The position of OZ at 
the moment ^ = is OZq. 
The motion of a point P whose abscissa is 

X = ne"^* cos (nt + a) 

where a, k, n, a are constants (a and k positive) is called a damped 
harmonic motion. It may be looked upon as a harmonic motion^ 
y^ose amplitude is decreasing. To study this motion let us 
again substitute a complex number 

z = 06""*' COS (7U+ a) + a^* sin (frf + a)t, 
or 

z = a6"**-6^'**^>% 
or 

where Zo is written for the complex constant ae^. 
The product 

is a complex number corresponding to a point Zi on the same 
radius as Zf^, coincident with Zq at the moment ^ = but apn 
proaching in a geometrical ratio after < = 0. In unit of time 
the distance of Zi from decreases in the constant ratio e"* : 1. 
The multiplication with 6"** turns OZi round through an angle 
rd. We may therefore describe the motion of Z as a uniform 
revolution of OZ round 0, Z at the same time approaching 
at a rate uniform in this sense that in equal times the distance 
is reduced in equal proportions (Fig. 26). At the moment ^ = 
the position coincides with Zq. We speak of a period of this 
motion meaning the time T = 2^/n in which OZ performs an 
entire revolution round 0, although it does not come back to its 
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original position. Any part of the spiral curve described by Z 
in a given time is geometrically similar to any other part of the 
curve described in an interval of equal duration. For suppose 
the second interval of time hap- 
pens r units of time later, we 
shall have for the first interval 




and for the second interval 

Now if z\ and %i are the values 
of z at two moments i\ and U of 
the first interval and z\ and Z/i 
the corresponding values of z' 
at the moments <i + t and fe + t of the second interval, we have 



'" t=8:p4 



Fig. 26. 



5 — ^*(li-li).^n(«i-««)» — 5l 
Zii Zi 



Therefore the triangle Z1OZ2 is geometrically similar to the 
triangle Z1OZ2. As Zi and Z2 may coincide with any points 
of the first part of the curve, the two parts are evidently geo- 
metrically similar. 

The projection of Z on the axis of x performs oscillations 
decreasing in amplitude. The turning-points correspond to those 
points of the spiral curve described by Z, where its tangent is 
parallel to the axis of y, that is to say, where the abscissa of dzjdt 
vanishes. 

Now 

dz 



or 



^ = zo(— i + n%)e^e'^ = (— i + ni)z 



dz 
dt 

z 
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where p and X are the modulus and angle of the complex number 

— A; + ni. 

Consequently, if we represent dz/dt by a point Z\ the triangle 

Z'OZ will remain geometrically similar to itself. The turning 

points of the damped oscillations correspond to the moments 

when OZ' is directed vertically upward or downward or when the 

angle of dz/dt is equal to t/2 or 3ir/2. The angle of z will then 

be x/2 — X or 3x/2 — X plus or minus any multiple of 2x. As 

the angle of z, on the other hand, is changing in time according 

to the formula 

tU+ a, 

we find the moments where the movement turns by the equation 
or 

frf + a = 3ir/2 - X + 2Nir, 

N denoting any positive or negative integral number. The time 
between two consecutive turnings is therefore equal to ir/n, that 
is, equal to half a period. All the points Z corresponding to 
turning points lie on the same straight line through the origin 
forming an angle 3ir/2 — X with the direction of the positive axis 
of X, The amplitudes of the consecutive oscillations therefore 
decrease in the same proportion as the modulus of z, that is 

to say, in half a period in the ratio ^"S"- 

Let us consider the vibrations of a system possessing one 
degree of freedom when the system is subjected to a force varying 
as a harmonic function of the time and let us limit our considera- 
tions to positions in the immediate neighborhood of a position 
of stable equilibrium. If the quantity x determines the position 
of the systeni^ the oscillations satisfy a differential equation of 
the form 

m ^ + fc-^ +n^x = Fcos (jrf)^ 
where m, k, n, p, F are positive constants. 

^ See for instance Rayleigh, Theory of Sound, Vol. I, chap. Ill, § 46. 
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This is another case where the introduction of a complex 
variable 

z = x + yi 



and the geometrical representation of complex numbers helps to 
form the solution and to survey the variety of phenomena that 
may be produced. 

In order to introduce z let us simultaneously consider the 
differential equation 

m^ + fe^ + w^y s= F sin p<, 

and let us multiply the second equation by i and add it to the 
first. We then have 

d Z , - dZ , n -n A* 

The movement of the point Z representing the complex number 
z then serves as well to show the movement corresponding to x. 
We need only consider the projection of Z on the axis of x. 

A solution of the differential equation may be obtained by 
writing 

Introducing this expression for z and cancelling the factor e^** 
we have 

2o(— WP^ + Jcpi + 71^) = Fy 

or 

F 

zo is a complex constant^ that may be represented geometrically 
as we shall see later on. 
This solution 

is not general. If 2' denotes any other solution so that 
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we find by subtracting the two equations 

^(z' -z) , ,d(z'- z) , ,, . . ^ 

or writing 

The general solution of this equation is 

where ui and ti2 are arbitrary constants and Xi and Xs are the 
roots of the equation for X 

mX^ + JfeX + n2 = 0, 



X2J 2m \4m2 



n^. 



If J^l^m ^ is greater than v?, so that the square root has a real 
value, Vk^/^m^ — n^ will certainly be smaller than k/2m. There- 
fore Xi and X2 will both be negative and the moduli of the complex 
numbers wie^** and n^e^ will in time become insignificant. If, 
on the other hand, F/4m^ is smaller than n^, both complex 
numbers uie^* and V/2fi^ correspond to points describing spirals that 
approach the origin, as we have seen above, in a constant ratio 
for equal intervals of time. Therefore they will also in time 
become insignificant. 
After a certain lapse of time the expression 

will therefore suffice to represent the solution. 

The point Z moves uniformly in a circle round of a radius 
equal to the modulus of zq, completing one revolution in the 
period 27r/p, the period of the force acting on the system. The 
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movement of the projection of Z on the axis of :e is given by 

a: = fo cos (p< + a), 

where fo is the modulus and a the angle of 20. It is a harmonic 
movement with the same period as that of the force? Fcospf, 
but with a certain difference of phase and a certain amplitude 
depending on the values of F, m, i, n, p. 

It is important to study this relation in order to survey the 
phenomena that may be produced. For this purpose the geo- 
metrical repre^ntation of complex numbers readily lends itself. 

In the expression for zo 

let us consider the denominator 

— tw]^ + Apt + n^ 

and let us suppose the period of the force acting on the system 
not determined, whfle the constants of the system m, A:, n and 
the amplitude of the force F have given values. The quantity p 
is the number of oscillations of the force during an interval of 
2tc units of time. This quantity p we suppose to be indeter- 
minate and we intend to show how the amplitude and phase 
of the forced vibrations compare with the amplitude and phase 
of the force for different values of p. 

Let us plot the curve of the points corresponding to the complex 
number 

where p assumes the values p = to + ^ • 

This curve is a parabola whose axi3 coincides with the axis of 
X and whose vertex is in the point a: = n^ y = 0. We find its 
equation by eliminating p from the equations 

« = n^ — mf^y y = ^> 
viz.. 
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But it is better not to eliminate p and to plot the different points 
for different values of p. In Fig. 27 the curve is drawn f or p = 

to 3 and the points for 
p^O,l, 2, 3 are marked. 
The ordinates increase in 
proportion to p; they are 
equal to 0, k, 2k, Sk for 
p = 0, 1, ?, 3. The dis- 
tance between the projec- 
tion of any point of the 
curve on the axis of x and 
It is equal to 0, m, 4m, 9m for 




Fio. 27. 



the vertex is proportional to p*. 

p = 0, 1, 2, 3.^ 

For any point P on the parabola let us denote the distance 

from by r and the angle between OP and the positive axis of 

ic by ^ so that 

n* — mp^ + kpi = re**. 
Then we have 



= *.^*< 



and consequently 



and 



zo = Te 



f 



(!>*-*)< 



a; = - COS (p< — ip). 

The amplitude F/r of the forced vibration is inversely propor- 
tional to r. Thus pur Fig. 27 shows us what the period of the 
force must be to make the forced vibrations as large as possible. 
It corresponds to the point on the parabola whose distance from 
is smallest. It is the point where a circle round touches the 
parabola. In Fig. 27 this point is marked R. It may be called 
the point of maximum resonance. When the constants of the 
system are such that the ordinate of the point, where the parabola 
intersects the axis of y is small in comparison with the abscissa 
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of the vertex, then OR wfll lie close to the axis of y (Fig. 28). In 
this case the angle between OR and the positive axis of x will be 
very nearly equal to 90**, that is to say, the forced oscillations will 
lag behind the force oscil- 
lations by a little less than a y 
quarter of a period. Keep- 
ing m and n constant, this 
will take place for small val- 
ues of k, i. e., for a small 
damping influence. A small ^iq. 28. 
deviation of p from the fre- 
quency of maximum resonance will throw the point P away from R, 
so that r increases considerably and <p becomes either very small 
(for values of p smaller than the frequency of maximum resonance) 
or nearly equal to 180** (for values of p larger than the frequency 
of maximum resonance). In other words for small values of k the 
maximum of resonance is very sharp. A deviation of the period 
of the force from the period of maximum resonance will lessen the 
amplitude of the forced vibration considerably. The lag of its 
phase behind that of the force will at the same time nearly vanish, 
when the frequency of theforce is decreased or itwill become nearly 
as large as half a period, when the frequency of the force is in- 
creased. For larger values of k the parabola opens out and this 
phenomenon becomes less marked. The minimum of the radius r 
becomes less pronounced. The angle between OR and the axis of 
z becomes smaller and smaller and for a certain value of k and all 
larger values the point R will coincide with the vertex of the para- 
bola. In this case, there is no resonance. When the period of 
the force increases indefinitely (p becoming smaller and smaller) 
the amplitude of the forced vibration will increase and will 
approach more and more to the limit 

F 

but there will be no definite period for which the forced vibra- 
tions are stronger than for all others. 



CHAPTER IL 



The Graphical Representation of Functions op One or 

More Independent Variables. 

§ 6. FuncHona cf One Independent Variable. — ^A function y of 
one variable x 

is usually represented geometrically by a curve, in such a way 
that the rectangular coordinates of its points measured in certain 
chosen units of length are equal to x and y. This graphical rep- 
resentation of a function is exceedingly valuable. But there is 
another way not less valuable for certain purposes, more used in 
applied than in theoretical mathematics, which here will occupy 
pur attention. 

Suppose the values of j/ are calculated for certain equidistant 
values of x, for instance: 

ar = - 6, - 5, - 4, - 3, -2, - 1, 0, 

+ l, + 2, + 3, + 4, + 5, + 6, 

and let us plot these values of ^ in a uniform scale on a straight 

line. Draw the uniform scale on 
one side of the straight line and 
mark the points that correspond 
to the calculated values of y on 
the other side of the straight line. 
Denote them by the numbers x 
that belong to them (Fig. 29). 
The drawing will then allow us to 
read off the value of y for any of 
p^Q^ 29. ^ *^^ values of x with a certain ac- 

curacy depending on the size of the 

scale and the nuinber of its partitions and naturally on the fine^ 
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ness of the drawing. It will also allow us to read off the value 
of y for a value of x between those that have been marked, if 
the intervals between two consecutive values of x are so small 
that the corresponding intervals of y are nearly equal. We can 
with a certain accuracy interpolate values of x by sight. On the 
other hand, we can also read off the values of x for any of the 
values of y. We shall call this the representation of a function 
by a scale. 

We can easily pass over to the representation of the same 
function by a curve. We need only draw lines perpendicular 
to the line carrying the scales through the points marked with 
the values of x and make their length measured in any given 
imit equal to the numbers x that correspond to them (Fig. 29). 

In a similar way we can pass 
from the representation of the 
function by a curve to the rep- 
resentation by a scale* 

The representation by a scale 
may be imagined to signify the 
movement of a point on a straight 
line, the values of x meaning the 
time and the points marked with 
these values being the positions 
of the moving point at the times 
marked. By passing over to the curve the niovement in the 
straight line is drawn out into a curve with the time as abscissa 
(Fig. 30). 

The representation by a scale is used in connection with the 
representation by a curve for the purpose of drawiijg a function 
of a function. ^ » r 

Let y be a function of x and x a function of t. Then we wish 
to represent y as a function of <. 

Let y — f(x) be given by a curve in the usual way^ and let 
X = <p{t) be given by a scale on the axis of x marking the points 
where / = 0, 1, 2, • • -,12. We then find the values of y corre- 
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spending to the values t = 0,l,2, • ■ -, 12 by drawing the ordi- 
nates of the curve y = f{x) for the abscissas marked ( = 0, 1, 2, 
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■ • ■, 12, These ordinates as a rule will not be equidistant. But 

as soon as we move them so as to make them equidistant, they 

form the ordinates of the curve 



t-» -:5 



:i i; 



-/WO) 



;i 



with t as abscissa (Fig. 31). 

The representation of a func- 
tion by a scale may be general- 
ised in the respect that neither 
of the two scales fadng one an- 
other on the straight line need 
necessarily be uniform. The in- 
tervals of both scales may vary 
^Q 02. from one side of the scale to the 

other. If the variation is suffi- 
dently slow the interpolation can nevertheless be effected with 
accuracy. We may look at this case as composed of two cases 
of the first kind. 

f{x) = y and y = gif). 
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These scales are placed together, so that the scale x touches the 
scale t 

while the scale y is cut out (Fig. 32). 

§ 7. The Principle of the Slide Ride. — ^Let us investigate how 
the relation between x and t changes by sliding the x- and ^-scales 
along one another. 

If we slide the a:-scale through an amount y = c so that a 

point of the or-scale that was opposite to a certain point y of the 

y-scale, now is opposite y + c, then the relation between x and t 

represented by the new position of the scales will be given by 

the equation 

m = g(t) + c. 

If X, t and x', If, denote two pairs of values that are placed 
opposite to one another, we shall have simultaneously 

m^g{t) + c, 

or by eliminating c 

m - g(t) = fix') - g(n. 

The ordinary slide rule carries two identical scales y = log a; and 
y == log t that are able to slide along one another, x and t running 
through the values 1 to 100. We therefore have 

log a: — log < = log x' — log f, 
or 

X ^ x' 

Fio. 33. 

That is to say, in any position of the x- and f-scale any two values 
X and t opposite each other have the same ratio (Fig. 33). This 
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is the principle on which the use of the slide rule is founded. 
It enables us to calculate any of the four quantities x, t, x\ t' 
if the other three are given. Suppose, for example, x^ t, xf 
known. We set the scales so that x appears opposite to i^ 









Fig. 34. 

then t is read oflF opposite to x'. On the other edges the slide 
rule carries two similar scales one double the size of the other 
(Fig. 34). We may write 

y = 2 log Z and y = 2 log T. 

By means of a little frame carrying a crossline and sliding over 
the instrument, we can bring the scales x and T or ^ and X op- 
posite each other. If, for example, for any position of the 
instrument x^ T and a;', T are two pairs of values opposite each 

other, then 

log a; - 2 log T ^logxf -2 log T, 
or 



X xf 



T2 



v2- 



If any three of the four quantities ar, T, x\ T are known the 
fourth may be read ofif. Thus we find the value 

xT^ 



by setting T opposite to a; and reading off the value opposite to 
T. Or we can find the value of 



^ 



X 



by setting X opposite to T and reading off the value opposite x\ 



GRAPHICAL BEPRESENTATION OF FXTNCTIONS. 45 

Let US reverse the part that carries the scales t, T so that x 
slides along T and X along t, but in the opposite order (Fig. 35). 




Fig. 35. 

The scales t, T may then be expressed by 

y = Z — log < and y = Z — 2 log T, 

I being the entire length of the scales. 

By setting the instrument to any position and considering the 

scales X and torX and T by means of the cross line we have 

loga:+ log« = loga:' + loe^' and logZ+ log T = logZ' + log F 

or 

xt = xr and XT = ZT, 

so that any two values opposite to one another have the same 
product. 

For X and T we have 

loga: + 21og r= log a;' + 2 log T, 
or 

xr = x'r\ 

Let us apply this to find the root of an equation of the form 

M^ + au = 6. 
Divide by w so that 

tr + a = - 
u 

and set r = 1 opposite to Z = 6. Then taking T = uwe find 
on the same cross line t^u^ and Z = bfu, so that we read the 
two values v? and b/u directly opposite to each other on the 
scales t and Z. If b/u is positive, it decreases while v? increases. 
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Running our eye along we have to find the place where the differ- 
ence hju — v^ is equal to a. Having found it the T-scale gives 
us the root of the equation. For example take 

w* — 5w = 3, 
or 

u 

We set r = 1 opposite X = 3 and run our eye along the scales 
X and t (Fig. 36), to find the place where ^ — 5 = -T. We find 






^Ki^l 



Of Oe M OT89Z0 9 ^ f f T 

Jg^ ► ! 1 ' '"lis i l-'S 8 



1 A i I tTTTI 



Fio. 36. 

it approximately at < = 6.2, and on the T-scale we read off 
T = 2.50 as the approximate value of the root. This is the 
only positive root. But for a negative root Zju is negative, 
and therefore the positive value of Zju plus v? would have to be 
equal to 5. We run our eye along and find t = 3.37 opposite to 
X = 1.63, approximately corresponding to T = 1.84. There- 
fore — 1.84 is another root. As the coeflScient of u^ in the first 
form of the equations vanishes it follows that the sum of the 
three roots must be equal to zero. This demands a second 
negative root approximately equal to — 0.66. To make sure 
that it is so, we set the instrument back and take the other end 
of the T-scale as representing the value T =• 1 and give it the 
position this end had before. Running our eye along the 
scales X and t, we find t = 0.43 opposite to X = 4.57, giving 
X+t= 5.00. On the T-scale we find 0.655, so that the third 
root is found equal to — 0.655. 

When b is negative there is always one and only one negative 
root. For u running through the values w = to — oo, v^—hju 
will run from — oo to + oo without turning. When h is positive 
there is always one and only one positive root; for then u^ — hju 
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• 

runs from — oo to + oo for i^ = to + oo. In the first case 

there may be two positive roots or none; in the second case there 

may be two negative roots or none. For positive values of a 

one root only exists in either case. This is easily seen in the first 

form of the equation 

w* + aw = 6, 

because from a positive value of a it follows that u^ + au will 
for w = — 00 to + 00, run from — oo to + oo without turning 
and will therefore pass any given value once only. 

In order to decide whether in the case of a negative value of a 
there are three roots or only one let us write 

1^* — - = — a. 
u 

For negative values of 6 we have to investigate whether there 
are positive roots. For positive values of u the function v^'-hju 
has a minimum, when the differential coefficient vanishes, t. e., for 

2i^ + -2 = 0, 
or 

u 

Having set our slide rule so that / gives us t? and X gives us 
— hjuy we find the value u where the minimum takes place by 
running our eye along and looking for the values X, t opposite 
each other for which X is twice the value of t 

2f = X. 

Then f + Z is the minimum of u^ — 6/w, so that there will be 
two or no positive roots according to ^ + X being smaller or 
larger than — a. For positive values of 6, we have to find out 
whether there are negative roots. The criterion is the same. 
After having set T = 1 opposite to h and having found the 
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positive root, we find the place where 

Then t + X is the minimum of all values that w* — b/u assumes 
for negative values of u. If the minimum is smaller than — a 
there are two negative roots; if it is larger there are none. If it 
is equal to — a the two negative roots coincide. 
For the equation 

ii»— 5 = -, 
u 

for instance, we find t = 1.31 opposite to Z = 2.62 (Fig. 36), 
so that 2t = 2.62 = X. Now t + X == 3.93 is smaller than 5, 
therefore u^ — 3/u will assume the value 5 for two negative 
values of u on either side of the value i^ = — T = — 1.143 
for which the minimum of u^ — 3/u takes place. 

On the same principle as the slide rule many other instruments 
may be constructed for various calculations. In all these cases 
we have for any position of the instrument 

fix) - git) = fix') - g(f), 

where x, t are any readings of the two scales opposite each other 
and x't' the readings at any other place, fix) and git) may be 
any functions of x and t It will only be desirable that they 
be limited to intervals of x and t, which contain no turning 
points. Else the same point of the scale corresponds to more 
than one value ot x or t and that will prevent a rapid reading 
of the instrument. 

Let us design an instrument for the calculation of the increase 
of capital at compound interest at a percentage from 2 per cent, 
upward. If x is the number of per cent, and t the number of 
years, the increase of capital at compound interest is in the pro- 
portion 

(' + m)'- 
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We can evidently build an instrument for which 



(i+ffi)'=(«+ffi)' 



For taking first the logarithm and then the logarithm of the 
logarithm, we obtain 

log< + loglogri + j^j = Iog<' + Ioglogfl+ ^j. 

We have only to make the ^c-scale 

y = + loglog(l + j^)- loglog(l + ~j, 

and the ^-scale 

y = log n — log t. 

For a; = 2 we have y = and therefore in the normal position 
of the instrument t = n. On the other end we have < = 1 and 
therefore y = log n. Now let us take n = 100, so that y = 2 
f or < = 1. Say the length of the instrument is to be about 24 
cm., then the unit of length for the y-scale would have to be 12 
cm. In the normal position of the instrument the readings x, t 
opposite to each other satisfy the equation 



(i+roo)*=(i+iio)" 



Opposite t = 1, we read the value Xi = 624 and this gives us 



( 



2 y^ xi 

l+ioo) =1 + 1^=1 + 6.24=7.24. 



A capital will increase in 100 years at two per cent, compound 
interest in the proportion 7.24 : 1, Or we may also say the 
number xi = 624 read ofiF opposite t = 1 is the amount which is 
added to a capital equal to 100 by double interest of 2 per cent, 
in 100 years. The same position of the instrument gives us the 
number of years that are wanted for the same increase of capital 
5 
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at a higher percentage. For all the values x, t opposite to each 
other satisfy the equation 



(l + i^) = 7-2^ 



For any other given percentage x and any other given number 
of years t the increase of capital is found by setting x opposite 
to t and reading the a:-scale opposite to ^ = 1. The only restric- 
tion is that the ratio is not greater than 7.24^ else ^ = 1 will 
lie beyond the end of the a:-scale. 

For a given increase of capital the instrument will enable us 
either to find the number of years if the percentage is given, or 
the percentage if the number of years is given, subject only to 
the restriction mentioned. 

We can build our instrument so as to include greater increases 
of capital by choosing a larger value of n. n = 1000, for in- 
stance, will make y = 3 for < = 1. If the instrument is not to 
be increased in size the scales would have to be reduced in the 
proportion 2 : 3. 

Let us consider another instance 

1 1 1 

y^x^ y^n^v 

In the normal position of the instrument the scale division 
marked x = oo corresponds to y = and is opposite to ^ = ti. 
If we have < = oo on the other end, the length of the instrument 
will correspond to y = 1/w. Let us choose n = 0.1, so that the 
length of the instrument is y = 10. That is to say, the unit of 
length of the y-scale is one tenth of the length of the instrument. 
For any position of the instrument we have 

1+1=1+1. 

If the scale division marked x = oo is opposite tot = c we can 
write a/ = oo, ^ = c and have 
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The instrument will therefore enable us to read off any one of 
the three quantities x, t, c, if the other two are given, the only 
restriction being that all three lie within the limits 0.1 to oo. 
The instrument may be used to determine the combined resistance 
of two parallel electrical re- 
sistances, for the resistances 
satisfy the equation 

R Ri R2 

Similarly it may be used 
to calculate the distances of an object and its image from the 
principal planes of any given system of lenses. For if / is the 
focal length and x and t the distances of the object and its im- 
age from the corresponding principal planes (Fig. 37), the equa- 
tion is 

On the back side of the movable part of an ordinary slide rule 
there generally is a scale 

y = 2 + log sin t 

When this part is turned round and the scale is brought into 

contact with the scale 

y = log X, 

we obtain for any position of the instrument 

log « — log sin ^ = log a:' — log sin /', 

x" 



or 



X 



sin t sin t 



f f 



for any two pairs of values x, t that are opposite each other. 
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Given two sides of a triangle and the angle opposite the larger 
of the two the instrument gives at once the angle opposite the 
other side. Similarly when two angles and one side are given, 
it gives the length of the other side. 

If x' = a is the value opposite to i' = 90**, we have 

a; = a sin i. 

Thus we can read the position of any harmonic motion for any 
value of the phase. 
An instrument carrying the scales 

y = log sin x and y = log sin i 

enables us to find any one of four angles ar, /, a:', H for which 

sin X sin a;' 



sin i sin i' 

if the other three are given. Thus, knowing the declination, 
hour angle and height of a celestial body, we can read the azimuth 
on the instrument. We have only to take x = 90® — height, 
i = hour angle, x' = 90** — declination, then f = azimuth or 
180® - azimuth. 

It is not necessary to carry out the subtraction 90® — height and 
90® — declination. The difference may be counted on the scale 
by imagining 0® written in the place of 90®, 10® in the place of 
80® and so on and counting the partitions of the scale backwards 
instead of forward. 

§ 8. Rectangtdar Coordinates with Intervals of Varying Size. — 
The two methods of representing the relation between two 
variables either by a curve connecting the coordinates or by 
scales facing each other lead to a combination of both. 

Suppose the rectangular coordinates x and y are functions of 

u and V, 

X = <p{u) and y = \l/(v). 

The function x = <p(u) is represented by a uniform scale for x 
on the axis of abscissae facing a non-imiform scale for u. The 
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function y = i^(v) is represented by a uniform scale for y oa 
the axis of ordinates fadng a non-uniform scale for v. Through 
^e scale-divisions u let us draw vertical lines, and through the 
scale-diviaons o let us draw horizontal lines. These two systems 
of parallel lines form a network of rectangular meshes of various 
sizes (Fig. 38), and any equation between u and o may be repre- 
sented by a curve in this plane. 

The usefulness of this method will be seen by some examples. 
It enables us by a dever choice of the functions y(u) and ^(e) 




to Amplify the form of the curve. It is ea^y seen, for instance, 
that a curve representing an equation /(w, o) = may always be 
replaced by a straight line, if we choose the u-scale property. 
For when the points u = I, 2, 3, 4, • • • of the curve are not on 
a straight line, let them be moved to a straight line without 
altering their ordinates (Fig. 39). This will change the «-scaIe 
but it will not alter the equation /(«, u) = now represented by 
the stnught line. 
Suppose we want to represent the relation 

o' ^ 6^ *■ 



where a and h are pven numbers. If « and t were ordinary 
rectangular coordinates the curve would be an ellipse. But if we 
make 

x = 1^ and y=^ 
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the equation of the line in rectangular coordinates becomes 

and the curve will therefore be a straight line running from a 
point on the positive axis of a; to a point on the positive axis of 
y. The point on the axis of x corresponds to the value u= ^ a 

on the ^^-scale, and the 
point on the axis of y cor- 
responds to the value v = 
sb 6 on the «-scale (Fig. 40). 
Any point in the straight 
line corresponds to four 
combinations +w, -^v^—Uy 
'{' v\ Uy — v\ — Uy ^ V, be- 
cause X has the same values 
for opposite values of u 




xfScCwj 



Fia.40. 



and y for opposite values of v. We can read « as a function of 
2^ or t^ as a function of v. 
If a second equation 



u 



2 



+ ri-=l 



bi 



2 



is given, we find the common solutions of the two equations by 
the intersection of the corresponding straight lines. Fig. 40 
shows the solutions of the two equations 



22+32-1 



and 



,2 



42 + 52 - h 



approximately equal to w = =*= 1.2 and v = ^ 2.4. 
Another function much used in mathematical physics 



tj = oe »« 
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may also be represented by a straight line by means of the same 
device. 
By making 

y = logv, X = u^, 
we obtain 

y=Iogo-^, 

where logv and log a are the natural logarithms of v and a. 
The li^scale is laid off on the axis of x and the i>-scale on the axis 
of y and we have to join the points u = 0, v = a and u = m, 

V = a/e. The point v = a/e is found by laying off the distance 

V = 1 to tJ = e from V =^ a downward (Pig. 41). We are not 
obliged to take the same units of length for x and y. 







0JiA5J2 A5±^ dJ!5 ^.«»" ^^ 



i>^a«'** 



Fig. 41. 



Suppose we had to find the constants a and m from two equa- 
tions 



and 






V2 = de mft' 



Our diagram would furnish two points corresponding to ui, vi 
and v^, V2. The straight line joining these two points intersects 
the axis of ordinates at i? = a and intersects the parallel through 
V = a/e to the axis of abscissae at t^ = m. 
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In applied mathematics the problem would as a rule present 
itself in such a form that more than two pairs of vfdues u, v 
would be given but all of them affected with errors of observation. 
The way to proceed would then be to plot the corresponding 
points and to draw a straight line through the points as best we 
can. A black thread stretched over the drawing may be used to 
advantage to find a strwght line passing as close to the points 
as possible (Fig. 42). 

In several other cases the variables u and e are connected with 
the rectangular co5rdinates x and y by the functions 
X = log u and y = log n. 




"Logarithmic paper" prepared with parallel lines for equidistant 
values of u and lines perpendicular to these for equidistant values 
of v is manufactured conmierciaUy (Fig. 43). 
By this device diagrams representing the relation 
uV = c, 
where r, *, c are constants are given by straight lines. For by 
taking the logarithm we obtain 

rx+ sy = log c. 

The straight line connects the point « = c''' on the u-8cale with 
the point » = (^ on the t>-scale. 

Logarithmic paper is further used to advantage in all those 
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cases where a variety of relations between the variables u and v 
are considered that differ only in u and v being changed in some 
constant proportion. If u and v were plotted as rectangular 
coordinates the curves representing the different relations be- 
tween u and V might all be generated from one of them by altering 
the scale of the abscissae and independently the scale of the ordi- 
nates, so that the appearance of all these curves would be very 

different. Let us write 

f{u, v) = 0, 

as the equation of one of the curves. The equations of all the 
rest may then be written 



(rs)-». 



where a, b are any positive constants. The points u, v of the 
first curve lead to the points on one of the other curves by taking 
w a times as great and 1^ 6 times as great. For if we write t^' = aw 
and v' = bu the equation f(u, r) = leads to the equation 
between m' and v': 



f{U)-o- 



Using logarithmic paper the diagram of all these curves be- 
comes very much simpler. The equation /(w, tj) = is equivalent 
to a certain equation <p(x, y) = 0, where a; = log w, y = log v. 
Now let x'y y' be the rectangular coordinates corresponding to 
w', ff so that 

x' = log 1^' = log w + log a = a; + log a, 
y' = log !^' = log tJ + log 6 = y + log h. 

The point «', y' is reached from the point x, y by advancing 

through a fixed distance log a in the direction of the axis of x 

and a fixed distance log h in the direction of the axis of y. The 

whole curve 

/(w, ^^) = 
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drawn on logarithmic paper is therefore identical with aU the 
curves 






It can be made to coincide with any one of the curves by 
moving it along the directions of x and y. 

§ 9. Functions of Two Independent Variables. — ^When a func- 
tion of one variable y = f(x) is represented by a ciu-ve, the values 
of X are laid off on the axis of x and the values of y are represented 
by lines perpendicular to the axis of x. In a similar way a 
function of two independent variables 

z = M y) 

may be represented by plotting x and y as rectangular coordinates 
and erecting lines perpendicular to the xy plane^ in all the 
points X, y, where /(a:, y) is defined and making the lengths of 
the perpendiculars proportional to z. In this way the function 
corresponds to a surface in space. Now there are practical 
difficulties in working with surfaces in space and therefore it 
appears desirable to use other methods, that enable us to represent 
functions of two independent variables on a plane. This may 
be done in the following way. 

Taking x, y as rectangular coordinates all the points for which 
f{Xy y) has the same value form a curve in the xy plane. Let 
us suppose a number of these feurves drawn and marked with the 
value of fix, y). If the different values of f{x, y) are chosen 
sufficiently close, so that the curves lie sufficiently close in the 
part of the xy plane that our drawing comprises, we are not only 
able to state the value of f{x, y) at any point on one of the drawn 
curves, but we are also able to interpolate with a certain degree 
of accuracy the value of f{x, y) at a point between two of the 
curves. As a rule it will be convenient to choose equidistant 
values of f{x, y) to facilitate the interpolation of the values 
between. The curves may be regarded as the perpendicular 
projection of certain curves on the surface in space, the inter- 
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sections of the surface by equidistant planes parallel to the 
xy plane. 

The method is the generaUzation of the scale-representation 
of a function of one variable. For a relation between i and x 
represented by a curve with t as ordinate and x as abscissa, is 
transformed into a scale .representation . by perpendicularly 
projecting certain points of the curve onto the axis of a, the 
intersections of the curve by equidistant lines parallel to the axis 
of X and marking them with the value of L A scale division in 
the case of a function of one variable corresponds to a curve in 
the case of a function of two independent variables. 

This method of representing a function of two independent 

variables by a plane drawing or we might also say of representing 

a surface in space by a plane drawing, is used by naval architects 

to render the form of a ship and by surveyors to render the form 

of the earth's surface and by engineers generally. Let us apply 

the method to a problem of pure mathematics. 

The equation 

z* + pz + g = 

defines z as a function of p and q. Let us represent this function 

by taking p and q as rectangular coordinates and drawing the 

lines for equidistant values of z. 

For any constant value of z we have a linear equation between 

the variables p and g, and therefore it is represented by a straight 

line. This line intersects the parallels p = 1 and p = — 1 at 

the points q=^ — :^ '- z and g = — ^ + z. Let us calculate 

these values for z = 0; =*= 0.1; =±: 0.2 • • • =±= L3 and in this way 

draw the lines corresponding to these values of z as far as they 

lie in a square comprising the values p = ^ 1 to + 1 and 

g = — 1 to + 1. Fig. 44 shows the result. On this diagram 

we can at once read the roots of any equation of the third degree 

of the form 

z» + pz + 3 = 0, 

where p and g lie within the limits — 1 to + 1. For p = 0.4 and 



60 



GBAPHICAL METHODS. 



g = — 0.2, for instance, we read z = 0.37, interpolating the value 
of z according to the position of the point between the lines 
z = 0.3 and z = 0.4. We also see that there is only one real 
root, for there is only one straight line passing through the point. 




Fig. 44. 

On the left side of the square there is a triangular-shaped region 
where the straight lines cross each other. To each point within 
this region corresponds an equation with three real roots. For 
example, at the point p=^ — 0.8 and g = + 0.2 we read z = 
— 1.00; + 0.28; + 0.72. On the border of this region two roots 
coincide. 

For values of p and q beyond the limits — 1 to + 1 the diagram 
may also be used. We only have to introduce z' = z/m instead 
of z and to choose m sufficiently large. 

Instead of 

2* + pz + g = 
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we obtain 



or dividing by m^ 



m^z^+prnzT + q^O, 



P ^, . 9 






m^ m* 



or 
where 



.^ 



z"^ + p'z' + g' = 0, 






By choosing a sufficiently large value of m, p' and q' can be 
made to lie within the limits — 1 to + 1 so that the roots z' 
may be read on the diagram. Multiplying them by m we 
obtain the roots z of the given equation. 

A function of two independent variables need not be expressed 
in an explicit form, but may be given in the form of an equa- 
tion between three variables 

g{u, V, w) = 0, 

and we may consider any two of them as independent and the 
third as a function of the two. The graphical representation 
may sometimes be greatly facUitated by modifying the method 
described before. The curves for constant values of one of the 
three variables, say Wy are not plotted by taking u and v as 
rectangular coordinates, but they are plotted after introducing 
new variables x and y, x b, function of u and y a function of v and 
making x and y the rectangular coordinates. 

In some cases, for instance, we can succeed by a right choice 
of the functions x = ip{u) and y = \l/{v) in getting straight lines 
for the curves w = const. This will evidently be the case, 
when the equation g(Uy i), w) = can be brought into the form 

a{w)(p(u) + b{w)\l/(v) + c{w) = 0, 

a, b, c being any functions of w, <p any function of u and ^ any 
function of v. 
For introducing 
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« = ^(w), y = ^W 

the equation will become 

oo: + fcy + c = 0, 

where a, b, c are constants for any constant value of w. 

As an example let us consider the relation between the true 
solar time^ the height of the sun over the horizon, and the declina- 
tion of the sun for a place of given latitude. Instead of the 
declination of the sun we might also substitute the time of the 
year, as the time of the year is determined by the declination of 
the sun. Our object then is to make a diagram for a place of 
given latitude from which for any time of the year and any 
height of the sun the true solar time may be read. 

In the spherical triangle formed by 
the zenith Z, the north pole P (if we sup- 
pose the place to be on the northern 
hemisphere) and the sun S (Fig. 45), the 
sides are the complements of the decli- 
nation 5, the height h, and the latitude 
ip. The angle t at the pole is the hour 
angle of the sun, which expressed in 
time gives true solar time. 
The equation between these four quantities may be written in 

the form 

sin A = sin ^ sin 5 + cos ^ cos S cos t 

The latitude <p is to be kept constant, so that t, Ji, S are the only 
variables. 
Now let us write 

X =^ cost, y = sin A, 

80 that the equation takes the form 

y = sin ^ sin 5 + x cos ^ cos 8. 

When X and y are plotted as rectangular coordinates, we obtain 
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a straight line for any value of 5. Let ns draw horizontal lines 
for equidistant values of A = to 90** and vertical lines for equi- 
distant values of < = — 180® to + 180** or expressed in time 
from midnight to midnight (Fig. 46). In order to draw the 




straight lines S = const.^ let us calculate where they intersect 
the vertical lines corresponding to a? = — 1 and a? = + 1 or 
expressed in time corresponding to midnight and to noon. For 
« = — 1 we have y = — cos (^ + 8), and for « = + 1 we have 
y .= cos (tp — 8). Let us draw a scale on the vertical a; = — 1 
showing the points y = — cos (^ + 8) for equidistant values of 
(^ + 8) and a scale on the vertical x ^ + 1, showing the points 
y = cos (^ — 8) for equidistant values of ^ — 8. The scale is 
the same as the scale for h, with the sole difference that the values 
of ^ — 8 are the complements of h and the values of ^ + 8 the 
complements of — h. For a latitude of 41**, for instance, we 
have 

For B ^ + 5^ — 5 

June21 23.6** 64.6* 17.5* 

September 23 and March 21 41* 41* 

December 21 -23.6* 17.6* 64.6* 
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The values »rf^ + 3 and p — 5 fumiah the interaectiona with 
the verticals a; = — 1 and x= -h i, so that the straight lines 
can be drawn corresponding to these daya of the year. The two 
outward lines are parallel but the middle hne is steeper. Their 
intersections with the horizontal line A = show the time of 
sunrise and sunset.* Strictly speaking the straight lines do 
not correspond to certain days. The straight line determined 
by any value of i changes its position continually as 5 changes 
continually. But the changea of 5 during one day are scarcely 
appreciable unless the drawing ia on a larger acale. 
If in the equation 

ax + by+c=0 

a and h are independent of w, only e being a function of w, all 
the straight lines w = const, are parallel. In thia case we are 
not obliged to draw the 
straight lines w = const. 
•* It will suffice to draw a 

line perpendicular to the 
lines w = const, and a 
scale on it that marks the 
points corresponding to 
equidistant values of w. 
On the drawing we place a 
00 sheet of transparent paper 
pja_ 47_ or celluloid,on which three 

straight lines are drawn is- 
suing from one point in the direction perpendicular to the u-acale, 
»-scaIe and i«-scale (Fig. 47). If we move the transparent material 
without turning it and make the first two lines intersect the u-and-o 
scale at given points, the wj-scale will be intersected at the point 
corresponding to the value of w. This method has the advantage 

I That b to Bay, the moment when the center of the sun would be Been on 
the horizon, if there were no atmospherical refraction. To take account of 
the refraction, the line ft — — 0.6° would have to be considered instead of 
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that we can use the same paper for a gteat many relations of 
three variables, as we can place a great many scales side by side. 
Or, in the case of one relation only, we may divide the region of 
the values u, v, w into a number of smaller regions and draw three 
scales for each of them, placing all the z^-scales or t>-scales or 
z^scales side by side. The drawing will then have the same 
accuracy as a drawing of very much larger size in which there 
is only one scale for each of the three variables. 

§ 10. Depiction of One Plane on Another Plane. — ^Let us now 
consider two quantities x and y each as a f imction of two other 

quantities u and v 

X = <p(u, v), 

y = ^{u, v). 

In order to give a geometrical meaning to this relation between 
two pairs of quantities let us consider x and y as rectangular 
coordinates of a point in a plane and u, v as rectangular coordi- 
nates of a point in another plane. We then have a corre- 
spondence between the two points. When the functions <p(u, v) 
and \l/(u, v) are defined for the values u, v of sl certain region, 
they will furnish for every point u, v of this region a point in 
the xy plane. Let us call this a depiction of the uv plane on 
the xy plane. Similarly a function of one variable x = <p(u) 
might be said to depict the u line on the x line. We may there- 
fore say that the depiction of one plane on another plane is, in 
a certain way, the generalization of the idea of a function of one 
variable. Let us suppose <p{u, v) and \l/(u, v) both to have only 
one value for given values of u and v for which they are defined. 
Then there will be only one point in the ocy plane corresponding 
to a given point in the uv plane. But to a given point in the 
Qcy plane there may very well correspond several points in the 
uv plane. 

Let us try to explain this by a graphical representation of the 
depiction of planes on each other. For this piupose we draw 
the curves x = const, and y = const, in the uv plane for equi- 
6 
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distant values of x and y. In the xy plane they correspond to 
equidistant lines parallel to the axis of x and to the axis of y. 
The point of intersection of two lines x = a and y = b corre- 
sponds to the points of intersection of the curves 
<fi(u, v) = a and i/(u, t) = b, 
in the ut plane. If in a certmn re^on of the uv plane, that 
we consider, they intersect only once there is only one point in 
the region of the uv plane considered and one point in the xy 
plane corresponding to each other. Jig. 48 shows the depiction 
of part of the uv plane on part of the xy plane. We have a net 
of square-shaped meshes in the xy plane and corresponding is a 
net of curvilinear meshes in the uv plane. 

Let us condder the curves x = const, in the uv plane as the 
perpendicular projections of curves of equ^ height on a surface 
extended over that part of the uv plane. From any point P 
of the surface corresponding to the values u, v we proceed an 



infinitely small distance, u changing to « + du, i> to v+ dv and 
xtox+ dx, where 



du = cos ad», dv = sin cds, 

where ds signifies the length of the infinitely small line from 
u, » to « + du, v+ dv in the uv plane and a the angle its direc- 
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tion forms with the positive axis of x. Let PN be a straight Kne 
whose projections on the u and v axis are equal to d<p/du and 
d<p/dv and let us write 

— = f cos X, T" = *• sm X, 

du dv 

r being the positive length of PN and X the angle between its 
direction and the positive axis of x. Then we have 

OX = -T-du+ -^00= rda COS (a ^ X), 

or 

dx 

^ = r cos (a - X). 

dx/ds measures the steepness of the ascent. It is positive when 
the direction leads upward and negative when it leads downward 
and its value is equal to the tangent of the angle of the ascent. 

From the equation 

dx 

^ = r cos (a - X) 

we see that the ascent is steepest for a = X, where dx/ds = r. 
The line PN in the u, ivplane shows the perpendicular projection 
of the line of steepest ascent on the surface x = <p{u, v) and the 
length of PN measured in the same unit of length in which u and 
V are measured is equal to the tangent of the angle of the ascent. 
Let us call the line PN the gradient of the fimction <p(u, v) at the 
point u, V. The direction of the gradient is perpendicular to the 
curve <p(u, v) = const, that passes through the point u, v; for in 
the direction of the curve we have 

dx 

— = 
da "' 

and therefore 

a - X = =t 90^. 

If PN' is the gradient of the function \p(u, v) at the point u, v, the 
angle between PN and PN' must be equal to the angle formed 
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by the curves x = const, and y = const, that intersect at the 
point u, V, or equal to its supplement according to the angle of 
intersection that we consider. 

Suppose the gradients PN and PN' do not vanish in any of 
the points in the region of the uv plane that we consider and 
that their length and direction vary as continuous functions of 
u and V. Let us further suppose that the gradient PN' (com- 
ponents: d\l//du, d\l//dv) is for the whole region on the left side 
of the gradient PN (components: dtp/du, d<p/dv), or else for the 
whole region on the right side of the gradient PN, then it fol- 
lows that any one of the curves x = const, and any one of the 
curves y = const, can only intersect once in the region considered. 

This may be shown by considering the directions of the curves 
X = const, and y = const, in the uv plane. Let us consider 
that direction on the curve y = const, in which x increases. If 
this direction deviates from PN the deviation must be less than 
90°, because dx/ds and therefore cos (a — X) is positive. Let us 
further consider that direction on the curve x = const, in which 
y increases. If it deviates from the direction of PN' the devia- 
tion must be less than 90°. Let us call these directions the 
direction of x (on the curve y = const.) and the direction of y 
(on the curve x = const.). Now if the gradient PN' is on the 
left of the gradient PN the y direction must also be on the left 
of PN (for if it were on the right of PN being perpendicular to 
PN it would form an obtuse angle with PN') and therefore it 
must be on the left of the x direction (for if it were on the right, 
PN' being perpendicular to the x direction would form an obtuse 
angle with the y direction, which we have seen to be impossible). 
Similarly it may be seen, that if PN' is on the right of PN, the 
direction of y will also be on the right of the direction of x. If 
therefore PN' is on the same side of PN in the whole region 
considered, the direction of y will also be on the same side of the 
direction of x for the whole region considered. This excludes 
the intersection of two curves x = const, and y = const, in more 
than one point. For, suppose there are two points of inter- 
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section and we pass along the curve y = const, in the direction of 
X. At the first point of intersection we pass over the curve 
X = const, from the side of smaller values of x to the side of 
larger values of x. Now if the values of x go on increasing 
as we go along the curve y = const, we evidently cannot get 
back to a curve x = const, corresponding to a smaller value of x. 
The only possibility of a second point of intersection would be 
that the direction in which the value of x increases on the curve 
y = const, becomes the opposite, so that in advancing in the 
same direction in which we came x would decrease again. 

The same holds for the curve 
a:= const. If we pass from one 
point of intersection with a 
curve y = const, along a curve 
X = const, to a second point 
of intersection with the same 
curve the only possibility is 
that the direction of y also be- 
comes opposite. This is ex- 
cluded as in contradiction with yiq, 49. 
the direction of y being on the 
same side of the direction of x throughout the whole region (Fig.49) 

It will be useful to look at it from another point of view. Let 
us consider a point A in the ud plane corresponding to the 
values u, v and let us increase u and « by infinitely small positive 
amounts du and do, so that we get four points ABCD, forming a 
rectangle corresponding to the coordinates. 

AiUyV) B iu + duyV\ C \UyV-]r dv) D : t^ + dw, i? + dv. 

In the xy plane these points are depicted in the points A^ 
By C, D, the intersections of two curves u and u+ du with two 
curves v and v + dv (Pig. 50). 

The projections of the line AB in the iCy plane on the axes of 
coordinates are obtained by calculating the changes of x and y 
for a constant value of v and a change du in the valu6 of u 
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Sip 3iL 

dxi = jj^du, dyi = -^du. 

Similarly the projections of AC are obtained by calculating the 
changes of x and y for a constant value of u and a change do in the 

value of t 

S tp 3^ 

dxt= T:dv9 dyt^-r-dv. 



dv 



dv 



Denoting the lengths of AB and AC by dsi and d!^ and the angles 
that the directions of AB and AC form with the direction of the 




Fig. 60. 



positive axis of x (the angles coimted in the usual way) by 71 
and 72 we have: 

dxi = (foi cos 71, dpi = dsi sin 71 



and 
or 

and 

We may call 



dx2 = dstcos 72, dy% = dsi sin 72, 



dtp 



dsi d^ 



dsi 



T- = COS 7i"T~> T~ = sm 7i-T^ 
du du du du 



dtp ds2 dyj/ . d82 



dv 



du y\duj "^VW 



the scale of depiction at ^1 in the direction AB and 
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-WW) 



dv 



the scale of depiction at A in the direction AC. It is here under- 
stood that the uv plane is the original^ which is depicted on the 
xy plane. If we take it the other way the scales of depiction 
in the directions AB and AC are the reciprocal values du/dsi 
and dvjds^. 
The area of the parallelogram ABCD in the acy plane is 

<foi<fes sin (7» - 71) = (1^^- ^^^^^'^'>- 

According to the way in which the angles 72 and 71 are defined 
sin (72 — 71) is positive, when the direction AC points to the left 
of the direction AB (assuming the positive axis of y to the left 
of the positive axis of x), and sin (72 — 71) is negative, when AC 
points to the right. Now dvdv is equal to the area of the rectangle 
ABCD in the uv plane. Therefore the value of 

dudv dv du 

is the ratio of the areas ABCD in the two planes and its positive 
or negative sign denotes the relative position of the directions 
AB and AC m the xy plane. We may call this ratio the scale 
of depiction of areas at the point A. 

du dv dv du 

is called the functional determinant of the functions <p{u, v) and 
4^iu,v). 

We have found the scale of depiction of lengths in the direc- 
tions AB and AC. Let us now try to find it in any direction 
whatever. From any point A in the uv plane, whose coordinates 
are u and v, we pass to a point D dose by whose coordinates are 
u + Au, V + Av. In the ocy plane we find the corresponding 
points A and D with coordinates (Fig. 51). 
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^. a: = ip(u, v) jj, x+ Ax^ ip(u + Aw, v + Av) 
• y s ^{u, v) 'y + Ay=tl^(u + Au, v + Av) 

We expand according to Taylor's theorem, and writing for 
shortness 



M 



we find 






> ^ = 



dv 



^^=s 



d4^ 
du 



^» = 



dv 



Ax = ^uAu + ^»Ai^ + terms of higher order. 
Ay = 4^uAu + ^vAv + terms of higher order. 



ik» 



^^Un 



T — J. — • 



^?j — rr 



jp 
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The length of AD and the angle of its direction we denote by 
Ar and a in the uv plane and by As and X in the xy plane* 
The limit of the ratio As/Ar, to which it tends, when D approaches 
A without changing the direction AD is the scale of depiction 
at the point A in the direction AD. 
Writing 

Au = Ar cos a, 

Av = Ar sin a, 
we obtain 

Aa; = (ipu cos a+ fp^ sin a)Ar + terms of higher order, 
Ay = (^„ cos a + ^^ sin a)Ar + terms of higher order. 

Dividing by Ar and letting Ar decrease indefinitely, we have in 
the limit 

dx 

^ = ^tt cos a + fp^ sin a. 
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dy 

"jt^ 4^u COS a + ^« sin a. 

For dx/dr and dyfdr we may also write ds/dr cos X, ds/dr sin X, 

"j~ cos A = ^u COS a-i" <pv sm a, 

■T"sin X = ^u cos a + ^» sin a. 

These equations show the scale of depiction ds/dr corresponding 
to the different directions X in the x, ^-plane and a in the u, v^ 
plane. 

By introducing complex numbers we can show the connection 
still better. 

Let us denote 

dx , dy . ds .. 

zi = ^tt + 4^uh 

2?2 = 9» + ^»t. 

Multiplying the second of the two equations by i and adding 
both they may be written as one equation in the complex form: 

2 = Zi cos a + 2^2 sin a. 

The modulus of z is the scale of depiction of the uv plane at the 
point A in the direction a. The angle of z gives the direction in 
the xy plane corresponding to the direction a. For a = we 
have z = Zi and for a = 90®, z = z^. 
Let us substitute 

cos tt = 2 * sin a 5=^ ^. 

and write 

zi + Zi/i . Zi — Zi/i 

so that the expression fbr z becomes 
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z= aer^+ be~^. 

This suggests a simple geometrical constructioa of the complex 
numbers z for different values of a. The term oe^ is represented 
by the points of a drcle described by turning the line that 
represents the complex number a round 
the origin through the angles a= • • • 2ir. 
The term be~*^ is represented by the 
points of a drcle described by turning 
the line that represents h round the ori- 
gin in the opposite direction through the 
angles a = ■ • • - 2t (Fig. 52). The 
addition of the two complex numbers 
Pia. 62. '^ *^*^ ^ "^ ^°'' *"y value of a is easily 

performed. The points corresponding 
to the complex numbers z describe an ellipse, whose two princi- 
pal axes bisect the angles between a and b. This is easily seen 
by writing 

a = rie<''^^", b = rje^'^f ''*. 

ao corresponds to the direction bisecting the angle between a 
and h and ai denotes half the angle between a and h (positive or 
negative according to the position of a and h). 

z = riet-^'+">' + rse<'^'— ", 
or 

= (fi + fs) cos (a — «i) + (ri — rs) sin (a — ai)i. 

Denoting the coordinates of the complex number ze~^ by J and ^ 
we have 

, . , = cos (a - on) and " = an (a — ai), 

and consequently the equation of an ellipse 

(f.+'i)'^(f.->i)' '■ 
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This ellipse turned round the origin through an angle equal to 
ao gives us the points corresponding to z. The principal axes 
are 2(ri + r2) and 2(ri — f2) (Fig. 53). The construction of 



z%^ .^ 



-2.1^-- 




X 



FiQ. 53. 

Fig. 53 is obvious. After plotting zi and z^ we find Ts^ji and 
— TS/iji by turning AZ2 through a right angle to the right and to 
the left. From these points lines are drawn to Zi. The bisection 
of these lines give a and 6. 

The figure shows that in case a and h have the same modulus, 
the triangle — Z^ji, Zi, Z^ji becomes equilateral and AZi is per- 
pendicular to the line joining — Ziji and ^2/1. In this case AZi 
and AZ2 would have the same or the opposite direction. But as 
Zi = ^u H" i^uif 2^2 = ^» + 4/vh this would mean that <Pu4^v — ^v^u 
= 0. 

The radii of the ellipse (Fig. 53) measured in the unit used 
give the different scales of depiction corresponding to the dif- 
ferent directions in the xy plane. We might also say the ellipse 
is the image in the xy plane of an infinitely small circle in the 
UD plane, magnified in the proportion of the infinitely small radius 
to 1, with its center in A. 

Z\ corresponds to a = and Z2to a = 90® and for a = to 90® 
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Z moves on the ellipse from Zi to Z^ thiough the shorter way. 
— Zi corresponds to o = 180" and — ^ to a = 270°. Now we 
have shown above that a positive value of the functional deter- 
minant ipu'i't — PdV'u means that Zi is on the positive side of Zi, 
so that in this case Z moves in the posidve sense (that is, in the 
direction from the posiUve axis of x to the positive axis of y) with 
increasing values of a. With a negative value Z moves in the 
opposite direction. 

Let us now suppose that the curves x = const, and y = const, in 
the uv plane intersect except on a certain curve where their direc- 




tions coincide in the way shown in Pig. 54. On this curve the 
functional determinant D = ipu^v — Vv'l'u must vanbh because 
the directions of the gradients coincide. Let us see what the 
depiction on the xy plane is like. 

Running along one of the curves y = const., say y = y\, 
toward the curve Z> = we intersect the curves x = x^, Xi, Xt 
until at the point A on the curve a; = ari we reach the curve Z> = 0. 
In the a^y plane the corresponding path is a parallel to the axis 
of a: at a distance y\ passing through x^, xa, xi and reaching a 
point A at xi. If we now proceed on the curve y = yi'm the 
uv plane beyond the curve Z) = 0, we again intersect the curves 
Xi, xt, etc., but in the inverse order. Thus the corresponding 
path in the xy plane does not pass beyond A, hut turns ba<^ 
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through the same points x^y yi\ Xzy yi, etc. The same holds for 
any of the other lines y = const. If we trace the line in the 
Qcy plane that corresponds to the points in the uv plane, where 
the curves x = const, and y = const, touch, we find the depiction 
of the uv plane only on one side of the curve in the xy plane. 
The other side has no corresponding points u, v. However to 
every point C on this side of the curve, there are two correspond- 
ing points C in the uv plane, one on either side of the curve 
Z) = 0. Imagine two sheets of paper laid on the xy plane; let 
them both be cut along the curve AB. Retain only the two 
pieces on this side of the curve and paste them together along 
the curve. The uv plane is in this way depicted on the paper 
in such a way that there is one point and one only on the paper 
corresponding to each point in 
the region of the uv plane con- \ 

sidered. The curve D = in --^jfj 
the uv plane corresponds to the 
rim where the two pieces of pa- 
per are pasted together. Any 
line straight or curved passing 
over the curve D = in the uv 

plane,corresponds to a line running from one of the sheets onto the 
other. It need not change its direction abruptly when it reaches 
the rim and passes onto the other sheet. For it may touch the 
rim in the direction of its tangent. This is actually the rule 
and the abrupt change of direction is the exception. Any line 
LAL (Fig. 55) in the uv plane, whose tangent as it crosses the 
curve D = at A does not coincide with the common tangent 
of the curves x = const, and y = const, will correspond to a line 
in the ocy plane, that does not change its direction abruptly 
when it touches the rim. 

This is best understood analytically. Let us consider corre- 
sponding directions at the points A in the uv plane and in the 
Q^ plane. We have seen above that corresponding directions 
(Fig. 56) are connected by the equations 
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«i 



^ 



^ 



Fig. 50. 

ds dx 
COS Xt; = x~ ^«* c^s a + ^v sin a, 

sinXv; = ^ = \^u cos a + ^» sin a. 
At the point A we have 

^tt^« — ^r^tt = 0. 

Assuming that the gradients at A do not vanish, so that we 

can write 

^u = r cos 7 , fpv = r sin 7, 

ypu = r' cos 7', ^r = r' sin 7', 

where r and r' are positive quantities, the equation fpu^v-- <pv4^u=0 
reduces to sin (7 — 7O = 0, that is, 7 = 7' or 7 = 7'+ 180**. 
It follows therefore that : 



cos X-j- = f cos (a — 7), 



ds 



sin X^ = r' cos (a — 7') = ^ t' cos (a — 7). 

Consequently for all directions a in the uo plane for which 
cos (a — 7) is not zero, we have 



tgX= ^-. 

T 
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That is to say, we have in the xy plane only one fixed direction 
X and the opppsite corresponding to all the different directions 
a except only a direction for which cos (a — 7) = 0. In the 
latter case, that is, when the direction a is perpendicular to the 
direction. 7 of the gradient, i. «., in the direction of th6 curves 
X = const, and y = const., we have 

da 
cos X -T- = 0, 

ds 

sm Xt"= 0. 

dr 

Therefore ds/dr = and X remains indeterminate. Any direction 
X for which tg X differs from + r'/r corresponds to a fixed direction 
a = Y + 90°ora = 7 - 90°, while (fo/dr = 0. 

As the curve D = is depicted on the rim of the two sheets 
of paper, all those lines that intersect the ciu^e 2) = in a 
direction different from the direction of the curves x = const, 
and y = const, are depicted in the ary plane as curves having 
their tangent at A in common with the rim. All lines in one of 
the sheets of paper that touch the rim at ^ in a direction differ- 
ent from that of the rim must be the depiction of lines in the uv 
plane that reach A in the direction of the lines x = const, and 
y = const. The scale of depiction is zero in the direction of the 
curves x = const, and y = const. In any other direction a 
we find it different from zero for: 



I = V^(,3 + /*) cos^ (a - y). 

It is a maximum in the direction a = 7 or 7 + 180** perpendicular 
to the curves x = const, and y = const. 

It may help to understand all these details if we discuss an 
example where the depiction of the uv plane on the a^ plane 
has a simple geometrical meaning, the planes being ground plan 
and elevation of a curved surface in space. The rim in the 
xy plane is the outline of the surface, the projection of those 
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points where the tangential plane b perpendicular to the plane 

of elevation. 
Suppose a cylinder of circular section cut in two half cylinders 

by a plane through its axis. Suppose one of the half cylinders 

in such a position that its axis 
forms an angle d with the 
ground plan, the plan of ele- 
vation being parallel to its 
axis, Fig. 67. Let us intro- 
duce rectangular coordinates 
u, V in the ground plan and 
rectangular coordinates x, y 
in the plan of elevation. A 
point P on the cylinder is de- 
fined by certain values u, v 
which define its ground plan 
and certain values x, y which 
define its elevation. It is 
easily seen from Fig. 57 that 
we have 




and 



X ^ u 



y = t* tg8 H 1 Va^— tJ^, 

^ * cos 



where a is the radius of the section. Now let us consider the 
elevation of the points P as a depiction of their ground plan. 
The functions ^(w, v) and ^(w, v) in this case are 

ip{u, v) = w, 

1 



and 



cos ' 



(Pu =1, ipv — 0; ^t. = tg 8, ^„ = — 



V 



cos 



hVa^ - v' ' 



(Pu4^9 — (Pv4^U = — 



V 



COS 



SVa^-i^' 



GRAPHICAL REPRESENTATION OF FUNCTIONS. 



81 




Fig. 58. 



The functional determinant vanishes tor v = on the line EF. 
The lines y = const, are the intersections of the cylinder with 
horizontal planes. In the plan of 
elevation they are straight hori- 
zontal lines; in the ground plan 
they are ellipses (Fig. 58). As we 
pass along one of these curves we 
cross the line EF in the ground 
plan but we only touch it in the 
plan of elevation, retracing the hori- 
zontal line back again. The lines 
X = const, are straight lines in both 
planes, but in space they corre- 
spond to ellipses. Again as we 
cross EF in the ground plan we 
only touch it in the plan of eleva- 
tion and retrace the vertical line down again. Any curve on 
the cylinder that crosses EF in a direction not perpendicular to 
the plan of elevation is projected in the plan of elevation with 
EF as its tangent. For the real tangent in space lying in the 
tangential plane of the cylinder can have no other projection, if 
not perpendicular to the plan of elevation. In this latter case 

the projection of the tangent is a point 
and the tangent of the elevation is deter- 
mined by the inclination of the oscillatory 
plane. 

There is a particular case to be consid- 
ered, when the curve 7) = in the uv plane 
coincides with one of the curves x = const. 
OT y = const. (Fig. 59), assuming the gra- 
dients of the functions <p(u, v) and 4/(u, v) 
not to vanish at the points of this curve. We have seen that at 
a point where D = the scale of depiction must vanish in the 
directions of the curve x = const, or y= const. Let the curve 
D = coincide with a line x = const., then it follows that the 




Fig. 69. 
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length of the depiction of this curve is zero and the depiction 
must be contracted in a point. For the length of the depiction 
of a curve z = const, is given by an integral 



/ 



da 
dr 



dr, 




where dr denotes an element of the curve and ds/dr the scale 
of depiction in the direction of the curve. As da/dr is zero all 
along the curve the integral must necessarily vanish. 

As an example let us con- 
sider 

z = uv, 

y "= V. 

The lines z == const, in the uv 
plane are equilateral hyper- 
bolas, the lines y = const, are 
parallels to the axis of u (Fig. 
60). Along the axis of u we 
have at the same time y = 0, 
ar = and Z)= « = 0. The 
whole axis of u is depicted in 

the point ar = 0, y = of the zy plane. 

Let us finally consider the case where the scale of depiction 

at any point is the same in all directions, though it need not be 

the same at different points. 
Writing as before 

_ dx dy . _ ds^ 
dr dr dr ' 

the connection between the scale of depiction ds/dr and the 
angles X, a determining corresponding directions in the ocy plane 
and in the uv plane is given by the equation 

z = 2i COS a + Zzsin a, 
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or 

where 

a = i(2i + ^fi)> b = i(ai — Zi/i). 

In the ease where the scale of depiction ds/dr, that is to say, the 

modulus of z, is independent of a, one of the constants a or 6 

must vanish, as we see at once from the construction of z (Fig. 
52). Let us consider the case 6 = 0, 

z = ae^ = :r ^» 

dr \ 

The complex number a may be written | a \ e^, where | a \ 
denotes the modulus of a and ao the angle. Both may vary 
from point to point, but at every point they have fixed values. 
Consequently we have 

-7- = I a I and X = a + ao» 

That is to say, from an angle a determining a direction in the 
uv plane, we find the angle X determining the corresponding 
direction in the xy plane by the addition of a fixed value ao. 
Any two directions a, a' will therefore form the same angle as 
the corresponding directions X, X' in the xy plane. The same is 
true when a = and z = be~^. The only difference is that in 
this latter case the direction of z rotates in the opposite sense 
with increasing values of a. 

Analytically depictions of this kind are represented by func- 
tions of complex numbers, 

x+ yi = f(u + m) or x+ yi ^ f{u — m). 

Assuming the function to possess a differential coefficient we have 

dx dy. 

Sx dy. 
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and therefore dther 

2i = V* or 2i = — V». 

Hence in the first case 

a = K^i + aiA") = zu b= Kzi - %/») = 

and in the second case 

a = 0, b = Zi. 

§ 11. 0/Aer Methods qf Representing Relations between Three 
VaruMes. — The depiction of one plane on another may be used 
to generalize the graphical representation of a function of two 
variables or a relation between three variables, as we prefer 
to say. 

As we have seen before, an equation 

g(x, y,z) = 

between three variables x, y, z can be represented by taking x 
and y as rectangular coordinates and plotting the curves z = 
const. (Fig. 61) for equidistant val- 
ues of z. Suppose now the xy plane 
to be depicted on another plane. 
The lines x = const.,y = const, and 
z = const, will be represented by 
three sets of curves. The fact that 
three values x, y, z satisfy the equa- 
tion g{x, y, z) = Q is shown geo- 
metrically by the intersection of 
fm. 61. th* three corresponding curves in 

one point. 
Another method for representing certain relations between 
three variables u, v, w consists in drawing three curves, each 
curve carrying a scale. The values of u, v, w are read each on 
one of the three scales. The relation between three values u v 
ie is represented geometrically by the condition that the corre- 
sponding points lie on a straight line (Fig. 62). This method is 
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far more convenient than the one using three sets of curves. It is 
less trouble to place a ruler over two points u, v of two curves 
and read the value w on the scale of the third than to find the 
intersection of two curves u = const, and v = const, among sets 
of others, pick out the curve w = const, that passes through the 




Fia. 62. 

same point and read the value of w corresponding to it. For we 
must consider that the curves corresponding to certain values 
of u and v are generally not drawn, but must be interpolated and 
so must the curve w = const. It is true that interpolations are 
necessary with both methods, but the interpolation on scales 
hke those in Fig. 62 is easily done. 

It must however be understood that while the three sets of 
curves form a perfectly general method for representing any rela- 
tion between three variables, the other method is restricted to cer- 
tain cases. In order to investigate this subject more fully we 
shall have to explain the use of line coordinates. 

When we apply rectangular coordinates x, y to define a certain 
point in a plane, we may say that x determines one of a set of 
straight lines (parallel to the axis of ordinates) and y determines 
one of another set of straight lines (parallel to the axis of abscissas) 
and the point is the intersection of the two (Fig. 63, 7). A 
similar method may be used to determine a certain straight line 
in a plane. Let x determine a point on a certain straight line, 
X being its distance from a fixed point A on the line measured 
in a certain unit and counted positive on one side and negative 
on the other. Let y define a point on another straight line 
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parallel to the first, y being its distance from a fixed point B on 
the line measured in the same way as x. The straight line 
passing through the two points is thus determined by the values 





Fia. es, 

X and y and for all possible values of x and y we obtain all the 
straight lines of the plane except those parallel to the lines on 
which x and y are measured. For simplicity we choose AB 
perpendicular to the two lines (Fig. 63, IF). Let us call x and y 
the line coordinates of the line connecting the two points x and 
y in Fig. 63, //, in the same way as x and y in Fig. 63, /, are 
called the point coordinates of the point where the two lines 
X and y intersect. 

A linear equation between point coordinates 

y = mx + II 

is the equation of a straight line. That is to say, all the points 
whose coordinates satisfy the equation lie on a certain straight 
line. If, on the other hand, we regard x and y as line coordinates 
we find the analogous theorem: all the straight lines whose 
line coordinates satisfy the equation 

y = mx + /x 

pass through a certain point. The equation is therefore called 
the equation of the point. 

In order to show this let us first draw the line a; = 0, y = /i 
(APO in Fig. 64). If now for any value of x we make AR = x 
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FiQ. 64. 



and PQ = mx, the point of intersection of RQ and AP must be 

independent of x, for 

PO mx 

AO X 

The ratio PO/AO determines the position of and as it is 
independent of x and the positions of A and P are also inde- 
pendent of X, the same, is true for 0. 
For negative values of m, PO and 
AO have opposite directions so that 
lies between A and P. 

For a given point 0, we can find 
the corresponding values of m and /i 
by joining with the points A and 
the point corresponding to a; = 1. 
If P and Q are the intersections of 
these lines with the line on which y 
is measured, we have JSP = /* and PQ = m. Any point in the 
plane thus leads to an equation 

except the points on the line on which x is measured. For 
m = the equation reduces to 

that is, the equation of a point on the line on which y is measured. 

Instead of y = mx + fi, we might also write x = m'y + m'* 
and go through similar considerations changing the parts of x 
and y. This form does not include the points on the line on 
which y is measured, but it does include the points on the line 
on which x is measured. For these we have m' = 0. 

The general equation of a point in line coordinates is given in 

the form 

ax+ by + c ^ 0, 

from which we may derive either of the first-mentioned forms 
dividing it by a or 6. 
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Dividing by e another convenient form is obtained, 



"* +-^=1, 



— c — c 



or writing 



— e — e 

- + ^ = 1, 

Xq determining the point of intersection of the line BO (Fig. 64) 
and the x-line, while yo determines the point of intersection of 
the line AO with the y-line. 

A curve may be given by an equation 

ai(u)x + bi(u)y + Ci(u) = 0, 

in which ai(u), 6i(w), Ci(u) are functions of a variable u. Any 
value of u furnishes the equation of a certain point and as u 
changes the point describes the curve. Let us suppose the curve 
drawn and a scale marked on it giving the values of u in certain 
intervals sufficiently close to interpolate the values of u be- 
tween them. Two other curves are in the same way given by 

the equations 

a2(v)x + b2(v)y + Ci{v) = 0, 

(h{w)x + hz{w)y + cz{w) = 0, 

and scales on these curves mark the values of v and w. 

Now we are enabled to formulate the condition which must be 
satisfied by the values u, v, w in order that the three corresponding 
points lie in one straight line. If x and y are the line coordinates 
of the line passing through the three points, x and y must satisfy 
all three equations simultaneously. 

Consequently the determinant of the three equations must 
vanish 

ai(b2Cz — 63Q2) + 02(63^1 — biCz) + (hQ>iC2 — hci) = 0, 
and, vice versa, if the equation between w, Vy w may be brought 



GRAPHICAL REPRESENTATION OP FUNCTIONS. 89 

into this form where ai, 61, ci are any functions of u, 02, 62, (h any 

functions of v and a^, b^, cz any functions of w, we can form the 

equations 

aix + biy + ci = 0, 

a^ + b2y + C2 = 0, 

(h^+biy+ cz = 0, 

and represent them graphically by curves carrying scales for 
u, Vy w. The relation between u, v, w is then equivalent to the 
condition that the corresponding points on the three curves lie 
on a straight line. But it must be remembered that only a 
restricted class of relations can be brought into the required form, 
so that the method cannot be applied to any given relation. 
The equation of a point 

ax + by + c = 

remains of the same form, when the units of length are changed 
for X and y. If x' denotes the number measuring the same length 
as the number x but in another unit, the two numbers must have a 
constant ratio equal to the inverse ratio of the two units. There- 
fore, by changing the units independently, we have 

x=\x\ y = fiy\ 

and the equation of the point may be written 

d\x' + bfxy' + c = 0, 
or 

a V + by + c = 0, 

where a' = \a and b' = fjtb. 

It is sometimes convenient to define the line coordinates in 
another way. Let { and 17 denote rectangular coordinates 
measured in the same unit, then the equation of a straight line 
can be written 

V == tg <p^ + rio, 

where <p is the angle between the line and the axis of { and 170, 
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the ordinate of the point of intersection with the axis of 17. 
Now let us call tg ^ and 170 the line coordinates of the straight 
line represented by the equation and let us denote them by x 
and y. Thus the values of x and y define a certain straight line 
and any straight line not parallel to the axis of ordinates may 
be defined in this manner. The condition that a straight line 
X, y passes through a point i, 17 is expressed by the equation 

or 

y = — faj + 1?. 

If we fix the values of x and y, all the values {, 17 that satisfy this 
equation represent the points of the straight line x, y and we 
therefore call it the equation of the straight line. If, on the 
other hand, we fix the values of { and 17, all the values x, y that 
satisfy the equation represent the straight lines that pass through 
the given point {, 17, and therefore we call it the equation of the 
point. 
The more general form 

aa: + 6y + c = 
can be reduced to 

a c 

It therefore represents the equation of the point, whose rec- 
tangular coordinates are f = a/6 and 17 = — c/6. The case 

where 6 = or 

aa;+ c = 

represents the equation of a point infinitely far away in the 
direction ^ or the opposite direction <p + 180°, tp being defined by 

c 

tg ip= x= — - . 

All the straight lines, whose coordinates x, y satisfy the equation 

aa: + c = 
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correspond to the same value of a; but to any value of y. That 
13 to say, they are all parallel and all the straight lines of this 
direction belong to them. 

Let us now discuss some of the applications of line co5rdinate3 
to the graphical representation of relations between three 
variables. 

The relation 



may be written in the form 

log M + log r = log w, 
or 

a: + y = log w, 
when 

a; = log u and y = log v. 

Let us plot X and y as line co- 
ordinates on two parallel lines CF^g. 
65), with scales for the values of u Pio 65 

and 0. The equadons x = log u 

and y = logv may be regarded as the equations of the points of 
these two scales. The equation 

x + y = \ogw 

for any value of w is the equation of a point. It can easily be 
constructed as the intersection of any lines x, y satisfying its 
equation. For instance, the line x = log w, y = and the line 
x= 0, y = hg w. The first line is found by connecting the 
scale division m = w of the ix-scale with the point B, the second 
by connecting the scale division v = w of the n-scale with the 
point A, If the units of x and y are taken of the same length, the 
point of intersection will lie in the middle between the two lines 
carrying the u and c scales on a line parallel to the two other lines 
and the w-scale will be half the size of the other two (Fig. 65). 
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log U + log e = log W 
expresses the condition that the three equations 

a: = log «, y =logv, x + y= log w 

are satisfied simultaneously by the same values of x and y, that 
is to say, that the three points on the u, t, iv scales corresponding 
to the values of u, v, w Ue on the same straight line x, y. 
The more general relation 

u^ = v>, 
where a and fi are any given values, can be treated in the same 
manner. Thus the pressure and volume of a gas undergoing 
adiabatic changes may be represented. In this case we have 

JTB* = w, 
where p denotes the pressure, v the volume and k and w con- 
stants. 

For a given gas k has a given value, but w depends on the 
quantity of the gas considered. 

We write 

a: = log p, y = log v. 
The relation then takes the form 
x+ ky = log w, 

and represents a point which may be con- 
structed by the intersection of any two 
straight lines x, y, whose coordinates sat- 
isfy the equation, for instance 



and 



1 = 0, 
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The first line connects the point B (Fig. 66) with the scale 
division p = w of the p scale and the second line connects the 
point A with the scale division of the v scale for which y = k log w. 
A perpendicular from the point of intersection on AB meets it in 
0' and as the ratio AO'/O'B is 

equal to the ratio of the seg- +j|. ^^.i^o 

ments on the p and v scales 
log w/k log w — 1/k it is inde- 
pendent of w. All the points '^^^ ++o«5 
corresponding to different val- 
ues of w lie on the same par- 
allel to the p and v scales and 
the w scale may be obtained 
by a central projection of the 
p scale on this parallel from 
the center B (Fig. 66). We 
might dispense with the con- 
struction of the w scale as 
long as the straight line for 
the w scale is drawn. For in 
using the diagram we gener- 
ally start with values po, Vq 
and want to find other values 
p, V, for which 
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Fig. 67. 



The straight line connecting the scale divisions p and v intersects 
the w scale at the same point as the straight line connecting the 
scale divisions po and Vq, so that we need not know the value of 
po^o*- It suffices to mark the point of intersection in order to 
find the value of p, when v is given or the value of v when p is 
given. 

Another example is furnished by the equation 



tc^ + xw -}- y = 0. 
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If we regard x and y as line coordinates any value of w determines 
the equation of a point. We plot the curve formed by these 
points with a scale on it indicating the corresponding values of w. 
Any values of x and y determine a straight line whose inter- 
sections with the w scale f lU'nish the roots of the equation. Each 
point of the w scale may be constructed by the intersection of 
two straight lines, whose coordinates x, y satisfy the equation, 
for instance 

X ^ 0, y ^ ^ V? and a: = — w, y = 0.^ 

In Fig. 67 the w scale is shown for the positive values w = to 

w = 2.5. 

In the same manner a diagram for the solution of the cubic 

equation 

117* + a^ + y = 0. 

or of any equation of the form 

v/" + xvf" + y = 
may be constructed. 

§ 12. Relations between Four Variables, — ^The method can be 
generalized for relations between four variables. 
Suppose four variables u, v, w, t are connected by the equation 

giu, V, w, t) = 0, 

and let us assume that for any particular value t = to the resulting 

relation between u, v, w can be given by a diagram of the form 

considered consisting of three curves carrying scales for u, v and 

w. Let us further suppose that for other values of t the scales 

for u and v remain the same, but the scale for w changes. Then 

we shall have a set of w scales corresponding to different values 

of t Connecting the points that correspond to the same value 

of w we obtain a network of curves t — const, and w = const. 

(Fig. 68). Any two values u, v furnish a straight line intersecting 

1 For small values of tr, this combination is not good because the angle of 
intersection is small. One might substitute x =^2, y — — v^ — 2w tor the 
first line. 
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the network of curves. The points of intersection correspond to 
values of t and w that satisfy the given relation. 
Any relation of the form 

*.(«)/((, w) + <l>{v)g{t, w) + h{t, «.) = 

may be represented in this way, ip{v) denoting any function of 
u, V'(^) iiQy function of v and fit, 
w), g{t, w), h{t, w) any functions 
of t and w. 

In tMs case we need only in- 
troduce the line coordinates z, 
y, writing 

X = v(u), y = ^(e). 

We then obtain a linear equation 
between x and y, 

fit, w)x + g(t, w)y + hit, w) = o, 

which for any given vidues of t f,o, 68. 

and w represents the equation of 

a point. For a given value of t and variable values of w we obtmn 
a curve t = const, carrying a scale for w and for a series of vidues 
of ( we obtain a set of curves t = const. Similariy for a g^ven 
value of 10 and variable values of ( the equation furnishes a curve 
w = const., carrying a scale for ( and a series of values of w 
furnishes a set of curves w = const. From any given values 
of u and v the line codrdinates x and y are calculated and the 
points where this straight line defined by x and y intersects 
the network of the curves ( = const, and w = const, furnish 
the values /, w that satisfy the relation together with the given 
values of « and t. The relation between the height, azimuth, 
declination of a celestial body and the latitude of the point of 
observation may serve as an example. Let A, a, 5 denote the 
height, azimuth and declination and ip the latitude. The angles 
t/2 — tp, x/2 — A, ir/2 — 6 are the three sides of a spherical 
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Fio. 60. 



triangle PZS (Fig. 69) formed by the pole P, the zenith Z and the 
celestial body S. The azimuth is defined as the supplement of 

z the angle PZS. 

The equation is 

sin 5 = sin ^ sin A — cos <p cos h cos a. 

We write 

X = cos a, y = sin 5, 

so that the equation becomes 

y = sin ^ sin A — a: cos <p cos h. 

We shall in this case use the second system of line coordinates 

where x is the slope of the line measured by the tangent of the 

angle formed with the axis of 

abscissas and y is the ordinate 

of the intersection with the axis 

of ordinates. If {, ri denote the 

rectangular coordinates of the 

point, the equation of the points 

takes the form 



1? = a;{ + y or y = ri — ix, 
so that in our case we have 
i = cos (p cos h, 17 = sin <p sin h. 



The curves <p = const, and h = 
const, can be drawn by means 
of these formulas. It is easily 
seen that they are ellipses and 
that the curves <p = const, are 
the same as the curves h = const. 
For a definite value of <p and a 
variable value of h we find 




FiQ. 70. 



+ 



cos* <f> ' sin* <(> 



= 1, 
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and for a definite value of h and a variable value of <p 

1^ •_« L •■•• 



cos^ h sin^ h 

Any of the ellipses intersects all the others and in this way they 
form a network. A point of intersection of the ellipse <p = ci 
and the ellipse h = C2 also corresponds to the values A = ci and 
<p = C2, as the ellipse ^ = ci is identical with the ellipse h = ci 
and <p = 02 identical with h = 02 (Fig. 70). The easiest way to 
find this network consists in drawing the straight lines 

( + 17 = cos (^ — A), 

and perpendicular to them the straight lines 

{ — 17 = cos (^ + A), 

for equidistant values of ^ + A and ^ — A. The ellipses run 

diagonally through the rectangular meshes formed by the two 

systems of straight lines. The scales for <p and A are written 

on the axis of coordinates, both scales being available for both 

variables. The scale for 5 is written on the axis of ordinates 

and is identical with the scale for t and A on this axis. For the 

ordinate corresponding to a given value 5 = c is sin c, and this is 

also the ordinate of the point where the ellipse ^ = c or A = c 

intersects the axis of ordinates. The scale for the azimuth cannot 

be laid down in exactly the same way as that for <p, A and d 

because cos a determines the slope of the straight line z, y. 

Let us draw a parallel to the axis of ordinates through the point 

^ = 1, ri = and mark a scale for the azimuth on it, making 

17 = cos a (Fig. 70). A line connecting the origin with any scale 

division of this scale has the slope of the line x = cos a,y = sin 5. 

To bring it into the position of the line x, y it must be moved 

parallel to itself, until its point of intersection with the axis of 

ordinates coincides with the scale division 5. This suggests 

another way of using the diagram. Let a pencil of rays be 

drawn from the origin to the scale divisions of the azimuth scale 

(Fig. 70), and let it be drawn on a sheet of transparent paper 
8 
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placed over the drawing of the ellipses. For any given value 
of 8 it is moved up or down as the case may be so that the center 
of the pencil coincides with the scale division 8. As long as the 
celestial body does not materially alter its declination the dia^ 
gram in this position will enable us to find any of the three 
values <p, h, a from the other two. 

As a second example let us consider the relation between the 
declination 6, the azimuth a, the hour angle ^ of a celestial body 
and the latitude <p of the point of observation. 

The relation is found by eliminating the height h from the 

equation 

sin 5 = sin ^ sin h — cos <p cos h cos a. 

For this purpose we express sin h and cos h by the other angles 

and substitute these expressions for sin h and cos h. 

We have 

cos h = cos B sin //sin a, 

sin A = sin ^ sin d + cos <p cos d cos t 

Substituting these values we find 

sin 6= sin^ <p sin 5+ sin <p cos <p cos 8 cos /—cos ^ cos B sin t ctg a, or 

cos^ ^ sin 5 = sin ^ cos <p cos 5 cos t — cos ^ cos 5 sin t ctg a. 

Dividing by cos^ <p cos 5 we finally obtain 

sin t 
tg a = tg ^ cos / — -— — ctg a. 

cos <p 

In order to represent this relation graphically we introduce line 

coordinates 

X = ctg a and y = tg 8 

and find 

. . sin/ 

« = tg tf) cos / X. 

^ ® cos ^ 

Let us use the second system of line coordinates. The rec- 
tangular coordinates {, iy of the point represented by the equation 
are found from it equal to: 



GHAPHICAL HEPEESENTATION OF FUNCTIONS. 99 

. wnt , , 

* COS V ' " 

The curves r = const, are ellipses. 
The curves ( = const, are hyperbolas, 

_g___£__j 

ain^ ( cos' t 

The ellipses and hyperbolas are confocal, the foci coindding 
with the points £ = =»= 1, ij = 0, so that the curves intersect at 
right angles. 

The scale for ip may be written on the axis of ordinates at the 
points where it intersects the ellipses. It is identical with the 
scale for S, the ordinate in both cases 
being the tangent of the angle with the 
only difference that 5 is negative on 
the negative part of the axis and ip is 
not. The scale for t may be written 
on one of the ellipses corresponding to ■ 
the largest value of tp that is to be taken 
account of. This ellipse forms the 
boundary of the diagram, ao that 
la^er values of <p are not represented. 
Corresponding to the azimuth we draw 
a pencil of rays on a sheet of trans- 
parent paper, which is laid on the draw- 
ing of the curves. The center of the 
pencil is placed on the scale division 5 j.,g^ 7i 

and the azimuth is equal to the angles 

that the rays form with the positive direction of the axis of or- 
dinates (Fig. 71). It suflBces to draw the curves and the rays 
only on one side of the a^s of ordinates. At the apex of the 
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hyperbolas the value of i changes abruptly. The line ^ = 6^ is 
meant to start from the focus ( == 1, 17 = 0. When the center of 
the pencil of rays is in the origin the rays form the asymptotic 
lines of the hyperbolas, a = 15° corresponding to ^ = 1*, a =30® 
to ^ = 2^ and so on. 



CHAPTER IIL 

The Graphical Methods of the Differential and 

Integral Calculus. 

§ 13. Graphical Integration. — We have shown how the ele- 
mentary mathematical operations of adding, subtracting, multi- 
plying and dividing and the inverse operation of finding the 
root of an equation can be carried out by graphical methods and 
how functions of one or more variables may be represented and 
handled. But the graphical methods would lack generality and 
would be of very limited use, if they were not applicable to the 
infinitesimal operations of differentiation and integration. In- 
deed it is here that they are found of the greatest value. In 
many cases, where the calculus is applied to problems of natural 
science or of engineering, the functions concerned are given in a 
graphical form. Their true analytical structure is not known 
and as a rule an approximation by analytical expressions is not 
easily calculated nor easily handled. In these cases it is of vital 
importance that the operations of the calculus can be performed, 
although the functions are only given graphically. 

Let us begin with integration, because it is easier than differ- 
entiation and of more general application. 

Suppose a function y = f(x) given by a curve whose ordinate is 
y and whose abscissa is x. The problem is to find a curve, whose 
ordinate Y is an integral of the function /(a;). 



= £f(x)dx. 



Let us assume the unit of length for the abscissas independent 
of the unit of length for the ordinates. The value of Y measures 
the area between the ordinates corresponding to a and x, the 
curve y = f(x) and the axis of x in units equal to the rectangle 
formed by the units of x and y. 
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In the simple case where f(x) is a constant the equation 
y ss f{^) s c is represented by a line parallel to the axis of x and 



F = 1 cdx =s c[x — a). 




Fig. 72. 



Y is the ordinate of a straight line intersecting the axis of x at 
the point ar = a. The constant c is the change of Y for an 

increase of x equal to 1. 
If P is the point on the 
axis of a; for a; = — 1 and 
Q the point where the line 
y = c intersects the axis 
of ordinates (Fig. 72) the 
desired line is parallel to 
PQ. It is constructed by 
drawing a parallel to PQ 
through the point a; = a on the axis of x (Fig. 72, where a = 0). 
When a given value ci is added, so that the equation becomes 

y = c {x— a) + Ci 

it amounts to the same as when the straight line is moved in the 
direction of the axis of ordinates through a distance Ci. For 
a; = a we then have Y = Ci, so that we obtain the line 

Y — c{x — a) + Ci, 

by drawing a parallel to PQ through the point x = a, y = Ci. 

In the second place let us assume that the line y = f{x) consists 
of a number of steps, that is to say, that the function has different 
constant values in a number of intervals a; = a:i to a^, a:^ to Xz, 
etc., while it changes its value abruptly at a^, ans, etc. The 
line presenting the integral 



= frndx 



does not change its ordinate abruptly. It consists of a con- 
tinuous broken line, whose corners have the abscissas a^, Xz, etc. 
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The directions of the different parts are found in the way just 
described by the pencil of rays from P to the points a, fi, y, etc. 
(Fig. 73), where the horizontal lines intersect the axis of ordinates. 




Fig. 73. 

To construct the broken line we draw a parallel to Pa through 
the point x = xi (in Fig. 73 Xi is equal to 0) as far as the vertical 
X = X2. Through the point of intersection with the vertical 
a: = ir^ we draw a parallel to Pp as far as the vertical x = xz. 
Through the point of intersection with the vertical a; = «« we 
draw a parallel to Py and so on. 

Finally let us consider the case of an arbitrary function y = f(x) 
represented by any curve. In order to find the curve 



= £f(x)dx 



we substitute for y = f(x) a function consisting of different 
constant values in different intervals and changing its value 
abruptly when x passes from one interval to the next, so that 
the line representing this function consists of a number of steps 
leading up or down according to the increase or decrease of /(a:). 
These steps are arranged in the following way. The horizontal 
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part AiAt of the first step (Fig. 73) starts from any point Ai 
of the given curve. The vertical part AtBi and the following 
horizontal part BiB% are then drawn in such a manner that B1B2 
intersects the curve and that the integral of the given function 
as far as the point of intersection Kb is equal to the integral of the 
stepping line as far as the same point. That is to say, the areas 
between the stepping line and the given curve on both sides of 
the vertical part A^Bi have to be equal. When Kb is fixed the 
right position of A2B1 may be found by eye estimate. The eye 
is rather sensitive for differences of small areas. Besides a shift 
of A2B1 to the right or to the left enlarges one area and diminishes 
the other 90 that even a slight deviation from the correct position 
makes itself felt. In the same way the next step B2C1C2 is 
drawn with its vertical part B2C1 in such a position that the 
areas on both sides are equal. The integral of the given curve 
as far as Ke will again have the same value as that of the stepping 
line as far as Kc. And so on for the other steps. The integral 
of the stepping line is constructed in the way shown. It is 
represented by a broken line beginning at the foot of the ordinate 
of Ai. The comers lie on the vertical parts of the steps or 
their prolongations. It is readily seen that the broken line con- 
sists of a series of tangents of the integral curve 



= £f(^)d^> 



and that their points of contact with the integral curve lie on 
the same verticals as the points Ai, Kb, Kc, etc. (In Fig. 73 these 
points are denoted 0, 2, 3, • • • .) That these points lie on the 
integral curve follows from the arrangement of the steps which 
make the integral of the given function at Kb, Kc, • • • equal to the 
integral of the stepping line. Now in the points Ai, Kb, Kc "• 
the ordinates of the given curve coincide with those of the 
stepping line. Hence both integral lines must for these abscissas 
have the same direction. 

1 In Fig. 73 the lower limit is 0. 
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Having constructed the broken line and marked the points 
2, 3, 4, • • • (Fig. 73), the integral curve is drawn with a curved 
ruler so as to touch the broken line in the points, 0, 2, 3, • • •. 
As the given curve does not change its ordinate abruptly the 
integral curve does not change its direction abruptly. The 
drawing shows how well the integral curve is determined by the 
broken line. There is practically no choice in drawing it any 
other way without violating the conditions. 

The ordinate of the integral curve is measured in the same 
unit as the ordinate of the given curve y = S{x). It may some- 
times be convenient to draw the ordinates of the integral curve 
in a scale different from that of the ordinates of the given curve. 
For instance the value of the integral may become so large that 
measured in the same unit the ordinates of the integral curve 
would pass the boundaries of the drawing board, or else they may 
be so small that their changes cannot be measured with sufficient 
accuracy. In the first case the scale is diminished, in the latter 
case it is enlarged. This is done by altering the position of the 
point P, the center of the pencil of rays that define the directions 
of the broken line. If P approaches the directions Pa, Pp, • • • 
become steeper to the same degree as if keeping P unchanged we 
had increased the ordinates of -4 1-^2, 5i-B2, • • • in the inverse pro- 
portion of the two distances PO. Hence by diminishing the 
distance PO the ordinates of the resulting broken line are enlarged 
in the inverse proportion. On the other hand, by increasing the 
distance PO the ordinates of the resulting broken line are di- 
minished in the inverse proportion of the distancies, because the 
change of the directions Pa, Pj8, • • • caused Jby^ longer distance 
PO is the same as if the ordinates of A1A2, B1B2, • • • were di- 
minished in the inverse proportion. The broken line constructed 
by means of the longer distance P'O will therefore be the same as 
if the ordinates of the stepping line were diminished. It therefore 
leads to an integral curve whose ordinates are diminished in the 
same proportion (Fig. 74). 

The graphical integration of 
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Y = £f{x)dx 

is not limited to values x > a. The method is just as well applic- 
able to the continuation of the integral curve for x < a. The 




steps have only to be drawn from right to left. The lower limit 
a determines the point where the integral curve intersects the 
axis of X. 

There is a method for the construction of the vertical parts 
of the steps, which may in some cases be useful, though as a rule 
we may dispense with it and fix their position by estimation. 
Suppose that A and B (Fig. 75) 
are two points where the curve is 
intersected by the horizontal parts 
of two consecutive steps and that 
the curve between A and B is a 
parabola whose axis is parallel to 
the axis of x. The position of the 
P,Q 75_ vertical part of the step between A 

ondB can be then found bya simple 
construction. Through the center C of the chord AB (Fig. 75) 
draw a parallel CD to the axis of x, D being the point of inter- 
section with the parabola. The vertical part EH of the step in- 
tersects CD in a point whose distance from C is twice the distance 
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from Z). That tUs is the right position of EH is shown as soon 
as we can prove that the area ADBGA is equal to the rectangle 
EHBG. The area ADGBA can be divided in two parts^ the tri- 
angle ABQ and the part ADBCA between the curve and the 
chord. The triangle is equal to the rectangle FIBG, while ADBCA 
is equal to two thirds of the parallelogram MNBA, and hence 
equal to the rectangle EHIF. Both together are therefore equal 
to the rectangle EHBG, and the two areas between the stepping 
line and the curve on both sides of EH are thus equal. 

If the curve between A and B is sup- 
posed to be a parabola with its axis par- 
allel to the axis of ordinates the con- 
struction has to be modified a little. 
Through the center C of the chord AB 
(Fig. 76) draw a vertical line CD as far 
as the parabola. On CD find the point 
K whose distance from C is double the 
distance from D and draw through it a 
parallel to the chord AB. This parallel 
intersects a horizontal line through C at a point L. Then EH 
must pass through i. This may be shown in the following way. 
The area between the parabola ADB and the chord AB is equal 
to two thirds of the parallelogram MNBA, MN being the tan- 
gent to the parabola at the point D. If Z)' is the point of inter- 
section of NN and the horizontal line through C, we have evi- 
dently 

CL = iCD\ 

Therefore the rectangle EHIF is equal to the area ADBA be- 
tween the parabola and the chord and EHBG is equal to ADGBA. 
Any part of a curve can be approximated by the arc of a 
parabola with suflScient accuracy if the part to be approximated 
is sufficiently small. When the direction of the curve is nowhere 
parallel to the axis of coordinates^ both kinds of parabolas may 
be used for approximation, those whose axes are parallel to the 
axis of X and those whose axes are parallel to the axis of y. But 
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X 



Fio. 77. 



when the direction in one of the points is horizontal (Fig. 76), 
we can only use those with vertical axes and when the direction 
in one of the points is vertical we can only use those with hori- 
zontal axes. Accordingly we have to use either of the two con- 
structions to find the position of the vertical part of the step. 
Do not draw your steps too small. For, although the difference 
between the broken line and the integral curve becomes smaller, 
the drawing is liable to an accumulation of small errors owing 

to the considerable number 
of corners of the broken 
line and little errors of 
drawing committed at the 
comers. Only practical ex- 
perience enables one to find 
the size best adapted to 
the method. 

Statical moments of areas 
may be found by a double 
graphical integration. Let us consider the area between the curve 
y = /(^) (Fig- 77), the axis of x and the ordinates corresponding 
to a; = and x. = {. The statical moment with respect to the 
vertical through a: — { is the integral of the products of each 
element ydx and its distance f — a: from the vertical 

^ = J (€ — x)ydx. 
Let us regard If as a function of { and differentiate it: 

-^ = ((f - x)y)^^ + J^ d{ ^^ "" ^^^^ 

= 0+ J yd^. 

That is to say, a graphical integration of the curve y = f(x) 
beginning at a: = furnishes the curve whose ordinate is 

dM 
di • 
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Hence a second integration of this latter curve. will furnish 
the curve ikf as a function of f . As if vanishes for { = the 
second integration must also begin at the abscissa x = 0. 




Fig. 78 shows an example. Each ordinate of the curve found 
by the second integration is the statical moment of the area on 
the left side of it with respect to the vertical through this same 
ordinate. The ordinate furthest to the right is the statical 
moment of the whole area with respect to the vertical on the 
right. The statical moment of the whole area with respect to a 
vertical line through any point Xi is the integral 



jf 



{xi — x)ydx. 



Considered as a function of xi its diflferential coefiBcient is 

That is to say, the diflferential coefficient is independent of Xu 
hence the statical moment is represented by a straight line. As 
its diflPerential coefficient is represented by a horizontal line 
through the last point on the right of the curve 



J ydx, 
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the direction of the straight line is found by drawing a line 
through P and through the point of intersection Q of the hori- 
zontal line and the axis of ordinates (Fig. 78). The position of 
the straight line is then determined by the condition that 



X 



(xi — x)ydx 



for xi= (is equal to the statical moment 



M{i) = J^ (£ - ^)ydx. 



We have therefore only to draw a parallel to PQ through the 
last point R of the curve for i/({) found by the second integration. 
The ordinates of this straight line for any abscissa Xi represent 
the values of 

{xi — x)ydx 



i 



measured in the unit of length of the ordinates. The point of 
intersection E with the axis of x determines the position of the 
vertical in regard to which the statical moment is zero^ that is to 
say, the vertical through the center of gravity. 
The moment of inertia of the area 



X 



ydx 



about the axis a; = ( is found in a similar way. It is exoressed 
by the integral 

r = J (i-xyydx. 

Considered as a function of ( we find by differentiation 

= + 2j ({ - x)ydx. 
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That is to say, the differential coefficiefit is equal to double the 
statical moment about the same axis. This holds for every value 
of f. Hence we obtain ^7 as a function of { by integrating 
the curve for 3f ({). For { = we have T = 0, so that the curve 
begins on the axis of a: at {= 0. 
The integral 

ydx 



£ 



is zero for x ^ a. The curve representing the integral has to 
intersect the axis of a; at a; = a (admitting values of a: > a and 
X < a), and it is there that we begin the construction of the 
broken line. If instead we begin it at the point x = a, y = c, 
the only difference is that the whole integral curve is shifted 
parallel to the axis of ordinates by an amount equal to c upwards 
if c is positive, downwards if it is negative. But the form of the 
curve remains the same. It is different when this curve is 
integrated a second time. For instead of 



we now integrate 



J ydx 
J ydx + c. 



The ordinate of the integral curve is therefore changed by an 

amount equal to c(x — a) and besides if the second integral curve 

is begun &.t x = a, y = Ci instead of a: = a, y = the change 

amounts to 

c{x — a) + ciy 

so that the difference between the ordinates of the new integral 
curve and the ordinates of the straight line 

y = cix — a) + ci 

is equal to the ordinates of the first integral curve (Fig. 79). 

This effect of adding a linear function to the ordinates of the 
integral curve is also attained by shifting the pole P upward or 
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downward. For it evidently comes to the same thing whether 
the curve to be integrated is shifted upward by the amount c or 
whether the point P is moved downward by the same amount, so 
that the relative position of P and the curve to be integrated 
is the same as before. Changing the ordinate of P by — c adds 



/ 

/x-* 



^ 



y 



/ 



.1.1 y <^_ x"^ n 



.^'^ 




fff{x)dxdx 



fffix)dxdx 



Fig. 79. 

c(x — a) to the ordinates of the integral curve. c(x — a) is the 
ordinate of a straight line parallel to the straight line from the 
new position of P to the origin. 

By this device of shifting the position of P upward or down- 
ward the integral curve may sometimes be kept within the 
boundaries of the drawing without any reduction of the scale of 
ordinates. A good rule is to choose the ordinate of P about 
equal to the mean ordinate of the curve to be integrated. The 
ordinates of the integral curve will then be nearly the same at 
both ends. The value of the integral 



J ydx 



is equal to the difference between the ordinates of the integral 
curve and the ordinates of a straight line parallel to PO through 
the point of the integral curve whose abscissa is a. 
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When the ordinate of P is accurately equal to the mean ordinate 
of the curve to be integrated for the interval a: = a to 6 the 
ordinates of the integral curve will be accurately the same at the 
two ends. But we do not know the mean ordinate before having 
integrated the curve. 

After having integrated we find the mean ordinate for the 
interval a; = a to 6 by drawing a straight line through P parallel 
to the chord AB of the integral curve, A and B belonging to the 
abscissas x=a and x=b. This 
line intersects the axis of ordi- 
nates at a point whose ordinate 
is the mean ordinate. 

Suppose a beam AB is sup- 
ported at both ends and loaded 
by a load distributed over the 

beam as indicated by Fig. 80. That is to say, the load on dx is 
measured by the area ydx. Let us integrate this curve graph- 
ically, beginning at the point A with P on the line AB. The 
final ordinate at B 




Fig. 80. 



X 



ydx 



gives tne whole load and is therefore equal to the sum of the two 
reactions at A and B that equilibrate the load. Integrating this 
curve again we obtain the curve whose ordinate is equal to 



Y being written for 



J Ydx, 

J ydx 



The ordinate of this curve at any point a; = { represents the 
statical moment of the load between the verticals x = a and 
a; = J about the axis a; = J. Its final ordinate BM, Fig. 81, is 
the moment of the whole load about the point B, and as the reac- 
tions equilibrate the load it must be equal to the moment of the 
9 
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reactions about the same point and therefore opposite to the 
moment of the reaction at A about B. If the reaction at A is 
denoted by Fa we therefore have 

FJih - a) = r Ydx. 
That b to say. Fa is equal to the mean ordinate of the curve 



- j ydx 



in the interval x ^ a to b. The mean ordinate is found by 
drawing a parallel to AM through P which intersects the vertical 
through A at the point F so that AF « /*«• As DB is equal to 




Fio. 81. 

the sum of the two reactions a horizontal line through F will 
divide BD into the two parts BG = Fa and GD = Fj^ 

Shifting the position of P to P' on the horizontal line FG 
and repeating the integration 

' Ydx, 



f 



we obtain a curve with equal ordinates at both ends. If we 
begin at A it must end in B. Its ordinates are equal to the 
difference between the ordinates of the chord AM and the curve 
AM (Fig, 81), and represent the moment about any point of 
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the beam of all the forces on one side of the point (load and 
reaction). 

The area of a closed curve may be f omid by integrating over the 
whole bomidary. Suppose x = a and a; = 6 to be the limits of 
the abscissas of the closed curve, the vertical ar = a touching the 
curve at A and the vertical ar = 6 at B (Fig. 82). By A and B 
the closed curve is cut in two, both parts connecting A and B, 
Let us denote the upper part by y = fi(x) and the lower part 
by y = f2{x). The whole area is then equal to the diflFerence 




FiQ. 82. 



or equal to 



rfi{x)dx— I f2(x)dx, 
t/a 



We begin the integral curve over the upper part at the vertical 
a: = a at a point JB, the ordinate of which is arbitrary, and draw 
the broken line as far as F on the vertical a; = 6 (Fig. 82). Then 
we integrate back again over the lower part, continuing the 
broken line from F to G. The line EG measured in the unit of 
length set down for the ordinates is equal to the area measured 
in units of area, this unit being a rectangle formed by PO and 
the unit of ordinates. That is to say, the area is equal to the 
area of a rectangle whose sides are PO and EG. 
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The method is not limited to the case drawn in Fig. 82, where 
the dosed curve intersects any vertical not more than twice. A 
more complicated case is shown in Fig. 83. But in all those cases 




where the object is not to find the integral curve but only to find 
the value of the last ordinate the method, cannot claim to be 
of much use, because it cannot compete with the planimeter. 



For the construction of the broken line we have drawn the 
steps in such a manner that the areas on both sides of the vertical 
part of a step between the curve and the stepping line are equal. 
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It would have also been admissible to construct the stepping 
line in such a way that the areas on both sides of the horizontal 
part of a step are equal (Fig. 84). Only the broken line would 
consist of a series of chords instead of a series of tangents of the 
integral curve. The points Ka, Kb, • • • , where the horizontal 
parts of the steps intersect the curve would determine the ab- 
scissas of the points of the integral curve, where its direction is 
parallel to the direction of the broken line. But this forms very 
little help for drawing the integral curve. That is the reason 
why the former method where the broken line consists of a series 
of tangents is to be preferred. However where the object is only 
to find the last ordinate of the integral curve the two methods 
are equivalent. 

§ 14. Graphical Differentiation. — ^The graphical differentiation 
of a function represented by a curve is not so satisfactory as the 
graphical integration because 



v=/W 




Fio. 85. 



the values of the differential 
coeflBcient are generally not 
very well defined by the curve. 
The operation consists in 
drawing tangents to the given 
curve and drawing parallels 
through P to the tangents 
(Fig. 85). The points of in- 
tersection of these parallels 
with the axis of ordinates fur- 
nish the ordinates of the curve representing the derivative. 
The abscissa to each ordinate coincides with the abscissa of the 
point of contact of the corresponding tangent. The principal 
difficulty is to draw the tangent correctly. As a rule it can be 
recommended to draw a tangent of a given direction and then 
mark its point of contact instead of trying to draw the tangent 
for a given point of contact. A method of finding the point of 
contact more accurately than by mere inspection consists in 
drawing a number of chords parallel to the tangent and to 
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bisect them. The points of bisection form a curve that inter- 
sects the given curve at the point of contact (Fig. 86). When a 
number of tangents are drawn, their points of contact marked 
and the points representing the differential coefficient constructed, 

the derivative curve has to be 
drawn through these points. 
This may be done more accur- 
ately by means of the stepping 
line. The horizontal parts of 
the steps pass through the 
points while the vertical parts 
lie in the same vertical as the 
point of intersection of two 
consecutive tangents. The derivative curve connects the points 
in such a way that the areas between it and the stepping line are 
equal on both sides of the vertical parts of each step. Thus 
the result of the graphical differentiation is exactly the same 
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Fig. 87. 

figure that we get by integration, only the operations are carried 
out in the inverse order. 

A change of the distance PO (Fig. 87) changes the ordinates 
of the derivative curve in the same proportion and for the same 
reason that it changes the ordinates of the integral curve when we 
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are integrating, but in the inverse ratio. Any change of the or- 
dinate of P only shifts the curve up or down by an equal amount, 
so that if we at the same time change the axis of x and draw it 
through the new position of P the ordinates of the curve will 
remain the same and will represent the jdifferential coefficient. 

When a function /(a:, y) of two variables is given by a diagram 
showing the curves f(x, y) = const, for equidistant values of 
f(x, y) the partial differential coefficients can be found at any 
point XQf yo by means of drawing curves whose ordinates represent 
f(x, yo) to the abscissa x orf{xo, y) to the abscissa y and applying 
the methods explained above. For this purpose a parallel is 
drawn to the axis of Xy for instance, through the point xq, yo 
and at the points where it intersects the curves f(x, y) = const, 
ordinates are erected representing the values of /(a;, yo) in any 
convenient scale. A smooth curve is then drawn though the 
points so found and the tangent of the curve at the point Xo 
furnishes the diflferential coefficient dfjdx for x = Xo,y = yo* 

The diflferential coefficients df/dx, df/dy are best represented 
graphically by a straight line starting from the point x, y to 
which the differential coefficients correspond, and of such length 
and direction that its orthogonal projections on the axis of x 
and y are equal to df/dx and df/dy. This line represents the 
gradient of the function f(x, y) at the point a;, y} It is normal 
to the curve f{Xy y) = const, that passes through the point or, y, 
its direction being the direction of steepest ascent. Its length 
measures the slope of the surface z = /(a;, y) in the direction of 
steepest ascent. This is shown by considering the slope in any 
other direction. Let us change x and y by 

r cos a, T sin a 

and consider the corresponding change 

Az = f{x + r cos a, y + r sin a) — f{x, y) 

of the function. By Taylor's theorem we can write it 
^ See Chap. II, § 10. 
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— r COS a + T— r sin a + terms of higher order in r , 
ax oy 

a is the direction from the point x, y to the new point a: + r cos a, 
y + r sin a and r is the distance of the two points. Dividing 
Az by r and letting r approach to zero we find 

hm — = T- cos a + — sm a. 
r dx dy 

This expression measures the slope of the surface z = f(Qcy) 
in the direction a. Now let us introduce the length I and the 
angle X of the gradient, and write 

f- = / cos X, T" = i sin X. 
da: ' dy 

Then we have 

— cos a + T- sm a = / cos (a — X). 
dx dy 

That is to say, the slope in any direction a is proportional to 
c6s (a — X), it is a maximum in the direction of the gradient 
(a = X) and zero in a direction perpendicular to it and negative 
in all directions that form an obtuse angle with it. When all 
three coordinates are measured in the same unit, the length of 
I measured in this unit is equal to the tangent of the angle of 
steepest ascent. Hence the length of the gradient varies with 
the unit of length. When the unit of length in which the values 
of f(xy) are plotted is kept unaltered, while we change the unit 
of length corresponding to the values x and y, the length of the 
gradient varies with the square of the unit of length. 

§ 15. Differential Equations of the First Order. — In the problem 
of solving a diflFerential equation of the first order 

by graphical methods the first question is how to represent 
the differential equation graphically. If x and y are meant to 
be the values of rectangular coordinates, the geometrical meaning 
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of the differential equation is that at every point x, y, where 
f{x, y) is defined, the equation prescribes a certain direction for 
the curve that satisfies it. Let us suppose curves drawn through 
all those points for which j{x, y) has certain constant values. 
Each curve then corresponds to a certain direction or the opposite 
direction. Let us distinguish the curves by different numbers or 
letters and let us draw a pendl of rays together with the curves 
and mark the rays with the same numbers or letters in such a way 
that each of them shows the direction corresponding to the 







curve marked with that particular number or letter (Kg. 88). 
Our drawing of course only comprises a certain region in which 
we propose to find the curves satisfying the differential equation. 
It may be that fixy) is defined beyond the boundaries of our 
drawing. Those regions have to be dealt with separately. 

The graphical representation of the differential equation in 
the region considered consists in the correspondence between 
the curves and the rays. It is important to observe that this 
representation is independent of the system of coordinates by 
means of which we have deduced the curves from the equation 
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We can now introduce any system of coordinates {, 17 and find 
from our drawing the equation 

that is to say, we can find the value of ^((, 17) at any point £, 17 
of our drawing. If, for instance, the unit of length is the same 
for i and 17 we draw a Une through the center of the pencil of rays 
in the direction of the positive axis of { and a line perpendicular to 
it at the distance 1 from the center. The segment on the second 
line between the first Une and the point of intersection with one 
of the rays measured in units of length and counted positive in 
the direction of positive 17 fiu'nishes the value of ^(f, 17) for all 
the points {, 17 corresponding to that particular ray. In this 
respect the graphical representation of a differential equation 
is superior to the analytical form, in which certain coordinates 
are used and the transformation to another system of coordinates 
requires a certain amount of calculation. 

Now let us try to find the curve through a given point P on 
the ciu've marked (a) (Fig. 88) that satisfies the differential equa- 
tion. We begin by drawing a series of tangents of a ciu've 
that is meant to be a first approximation. Through P we draw 
a parallel to the ray (a) as far as the point Q somewhere in the 
middle between the curves (a) and (6). Through Q we draw a 
parallel to the ray (6) as far as R somewhere in the middle 
between the curves (b) and (c). Through R we again draw a 
parallel to the ray (c) and so on. The curve touching this 
broken line at the points of intersection with the ciu'ves (a), 
(6), ••• is a first approximation. But we need not draw this 
curve. In order to find a better approximation we introduce a 
rectangular system of coordinates x, y, laying the axis of x some- 
what in the mean direction of the broken line. Let us denote 
by yi the function of x that corresponds to the curve forming the 
first approximation. The second approximation jfe is then ob- 
tained as an integral curve of /(a:, yi), that is, of dyi/dx 
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denoting by x^p^ yp, the coordinates of P. For this purpose the 
curve whose ordinates are equal to /(a;, y\) or dyijdx has to be con- 
structed first. The values of fix, yi) are found immediately at 
the points where the first approximation intersects the curve 
(a), (6) • • • by differentiation in the way described above. A 
line is drawn through the center of the pencil of rays parallel to 
the axis of x and a line perpendicular to it at a convenient dis- 
tance from the center. This distance is chosen as the unit of 
length. The points of intersection of this line with the rays de- 
termine segments whose lengths are equal to the values o!f(x, yi) 
on the corresponding curves. These values are plotted as ordi- 
nates to the abscissas of the points where the first approximation 
intersects the curves (a), (6), • • • and a curve 

is drawn (Fig. 88). This curve is integrated graphicaUy begin- 
ning at the point P and the integral ciu've is a second approxi- 
mation. Again we need not draw the curve. The broken line 
suflSces, if we intend to construct a third approximation. In 
this case we have to repeat the foregoing operation. This can 
now be performed much quicker than in the first case because the 
values of f(x, y) on the curves (a), (6), • • • have already been 
constructed and are at our disposal. In order to find the curve 

Y = S{x, y^) 

we have only to shift the same ordinates to new abscissas and 
make these coincide with the abscissas of the points where the 
second approximation intersects the curves (a), (6), •••• The 

curve 

Y = fix, 2fe) 

is then drawn and integrated graphically, beginning at the point 
P. 
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Suppose now the integral curve did not differ from the second 
approximation, it would mean that 

»i = yp + I /(a^» y2)^f 

or that 

that is to say, that y^ satisfies the differential equation. 

If there is a perceptible difference the integral ciu^e represents 
a third approximation. It has been shown by Picard that pro- 
ceeding in this way we find the approximations (under a certain 
condition to be discussed presently) converging to the true solu- 
tion of the differential equation, so that after a certain number 
of operations the error of the approximation must become 
imperceptible. 

Denoting by yn the function of the nth approximation we have 

yn+i = yp + I /(^> yn)da:. 

The true solution with the same initial conditions y — yp for 
X =^ Xp satisfies the equation 

I f(x,y)cLc. 

Hence 



or 



r /(^> yn) - / (a?> y) / V. 



^ yn-y 

Let us now suppose that the absolute value of 

SJXy yn) - S{x, y) 

yn-y 

for all the values of x, y, yn within the considered region does 
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not surpass a certain limit M, then it follows that a certain relation 
must exist between the maximum error of yn, which we denote by 
Cn and the maximum error of yn+i, which we denote by Cn+i. 
The absolute value of the integral not being larger than 



Men \X — Xn 

(,| a; — a:n I denoting the absolute value of a; — ar„) we have 

CnJfl ^ M\X— Xn\ en. 

Hence as long as the distance x^ Xn over which the integration 
is performed is so small that 

M\x- Xn\ ^ A; < 1, 

h being a constant smaller than one, the error of y^^i cannot be 
larger than a certain fraction of the maximum error of yn. 
But in the same way it follows that the error of yn cannot be 
larger than the same fraction of the maximum error of yn-i, and 

so on, so that 

But as ei is a constant and k a constant smaller than one, k^ei 
must be as small as we please for a sufficient large value of n. 
That is to say, the approximations converge to the true solution. 
M being a given constant the condition of convergence 

M\x- xA <h<l 



limits the extent of oiu* integration in the direction of the axis of x. 
But it does not limit our progress. From any point P' that we 
have reached with sufficient accuracy we can make a fresh start, 
choosing a new axis of x suited to the new situation. As a 
rule it does not pay to trouble about the value of M and to try 
to find the extent of the convergence by the help of this value. 
The actual construction of the approximations will show clearly 
enough how far to extend the integration. As far as two consecu- 
tive approximations show no difference they represent the true 
curve. 
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Suppose that 

M Vn) - M y) 
Vn-y 

has the same sign for all values z^ y, yn concerned. Say it is 
negative. Suppose further that yn — y is of the same sign for 
the whole extent of the integration 

•/wp yn y 

that is to say, the approximative curve yn is all on one side of the 
true curve. Then if a? — aip is positive, yn-n — y must evidently 
be of the opposite sign from yn — y, or the approximative curve 
yn+i is all on the other side of the true curve from yn. For these 
and all following approximations the true curve must lie between 
two consecutive approximations. If the first approximation yi is 
all on one side of the true curve the theorem holds for any two 
consecutive approximations. This is very convenient for the esti- 
mation of the error. 

In Fig. 88 

/fe yn) - f(x, y) 

yn-y 

is negative from the point P as far as somewhere near S. The 
first approximation is all on the upper side of the true curve. 
Therefore the second approximation must be below the true 
curve at least as far as somewhere near S. 

When the sign is positive the same theorem holds for negative 
values of a; — Xp. If the integration has been performed in the 
positive direction of x, it may be a good plan to check the result 
by integrating backwards, starting from a point that has been 
reached and to try if the curve gets back to the first starting 
point. In this direction we profit from the advantage of the 
true curve lying between consecutive approximations and are 
better able to estimate the accuracy of our drawing. 

We have seen that the convergence depends on the maximum 
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absolute value of 

/(a?, Vn) - Six, y) 

Vn-y 

for all values of x, y, yn concerned. In order to find the maximum 
value we may as well consider 

for all values of a;, y within the region considered. For if we 
assume dfjdy to be a continuous function of y, it follows that 
the quotient of differences 

/fe Vn) - Six, y) 
Vn—y 

must be equal to dSldy taken for the same value of x and a value 
of y between y and yn. This is immediately seen by plotting 
Six, y) as ordinate to the abscissa y tor a fixed value of x. The 
value of the quotient of differences is determined by the slope 
of the chord between the two points of abscissas y and y„. The 
slope of the chord is equal to the slope of the curve at a certain 
point between the ends of the chord. The value of df/dy at this 
point is equal to the value of 

Si^y Vn) - Si^f y) 

yn-y 

Now let us consider how the coordinate system may be chosen 
in order to make dSJdy as small as possible and thus obtain the 
best convergence. For this purpose let us investigate how the 
value of dSJdy changes at a certain point, when the system, of 
coordinates is changed. 

Let us start with a given system of rectangular coordinates {, 
17 with which the differential equation is written 

The direction of the curve satisfying the differential equation 
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forms a certain angle a with the positive axis of { determined by 

tg a = ^ = <p(i, v) 

(assuming the coordinates to be measured in the same unit). 
Now let us introduce a new system of rectangular coordinates 
X, y connected with the system (, 17 by the equations 

« = { cos 0) + ^? sin «, 

y = — { sin « + 17 cos 0), 

which are equivalent to 

{ = ar cos M — y sin «, 

17 = a: sin « + y cos w, 

0) being the angle between the positive direction of x and the 
positive direction of {, counted from { towards x in the usual way. 
The angle formed by the direction of the curve with the positive 
direction of the axis of a; is a — «, and therefore 

^= tg (a- w) =f(x,.y). 

Consequently we obtain for a given value of o) 

^ 1 da 

dy cos^ (a — «) dy ' 

or remembering that a is given as a function of ( and 17, 

^ 

dy cos^ (< 

For simplicity's sake we shall assume that the axis of J is the 

tangent of the curve <p{i, 1?) = const, that passes through the 

given pointy so that da/d^ = 0. 

We then have 

df 1 da 

dy cos^* (a — «) di; 
and our object is to find how df/dy varies for different values of 



1 ( da . .da \ 

^ \ 'I — T^ sm 0) + — cos 6) I. 

,a — «) \ dj dri J 
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«. The value of dajdii is independent of «; it denotes the value 
of the gradient of a, which we represent by a straight line drawn 
from the origin A (Fig. 89) perpendicular to the curve a = const, 
or ^({, 7)) = const. 

It is no restriction to assume the value of da/drj positive; it 
only means that the direction of the positive axis of ly is chosen 




tf-6 



Fig. 89. 



in the direction of the gradient. Let us draw the line AB (Fig. 
89) in the direction of the positive axis of { and of the same length 
as the gradient. 

In order to show the values of df/dy for the different positions 
of the axis of x let us lay off the value of df/dy as an abscissa. 
For instance for <a = a, df/dy assumes the value 



da 

drj 



COS a. 



The abscissa corresponding to this value is AB' (Pig. 89), the 



10 
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orthogonal projection of AB on the axis of x. For any other 
position AC (Fig. 89) corresponding to some other value of u, 
we find dajdii cos a by orthogonal projection of AB on AC. Then 
the diviaon by cos (a — w) furnishes AC and a second division 
by cos (or — a) leads to AC. Hus a certwn curve can be 
constructed whose polar coordinates are r » tffldy and a, the 
equation in polar coSrdinates bong 

"" — ' -- da 

or [r cos (o — «)r= t- r cos w. 



di7 co^ (a — m) 
In rectangular coordinates {, q the equation assumes the form 



This shows that the equation is a parabola, the axis of which is 
perpendicular to the direction a. AB' is a chord and the gradient 




AG is a tangent of the parabola. Bisecting AB' in E, drawing EK 
perpendicular to ^£' as far as the axis of i; and bisecting EK in 
D, we find D the apex of the parabola. The three points A, B', D 
together with the gradient will suffice to give us an idea of the 
size and sign of Sfl^y for the different positions of the positive 
axis of X. 
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dfldy vanishes when the axis of x is perpendicular to the curve 
a » const.; so that it seems as if this were the most favorable 
position. We must^ however, bear in mind that the axis of x 
is kept unaltered for a certain interval of integration. When we 
pass on to other points the axis of x is no longer perpendicular 
to the curve a = const, there. The position of the axis of x is 
good when the average value of dfjdy is small. In Fig. 90 the 
parabolas are constructed for a number of points on the first 
approximation of a curve satisfying the differential equation. 

If we want to make use of the parabolas to give us the numerical 
values of dfjdy the unit of length must also be marked in which 
the coordinates are measured. The numerical value of dfjdy 
varies as the unit of length and therefore the length of the line 
representing it must vary as the square of the unit of length. 
But if we draw a line whose length measured in the same unit is 

equal to tttt > this line would be independent of the unit of 
length. For if Ms the line representing the unit of length and 
/', Z" the lines representing the values dfjdy and tttt- > ^jdy 

would be the ratio Vjl and ^ the ratio V'jl) hence V = P/Z'. 

Since /' varies as P with the change of the unit of length V is 
independent of the unit of length. This line Z" represents the 
limit beyond which the product 

becomes greater than 1. If dfjdy remained the same this would 
mean the limit beyond which the convergence of the process of 

approximation ceases. We might lay off the length of ^/ in 

the different directions in the same way as dfjdy has been laid 
off. The result is a curve corresponding, point by point, to the 
parabola, the image of the parabola according to the relation of 
reciprocal radii. But all these preparations as a rule would not 



132 GRAPHICAL METHODS. 

pay. It is better to attack the integration at once with an axis 
of X somewhat perpendicular to the curves a = const, as long 
as the direction of the curve forms a considerable angle with 
the curve a = const, and to lose no time in troubling about the 
very best position. The convergence will show itself, when the 
operations are carried out. When the angle between the direction 
of the curve that satisfies the differential equation and the curve 
a = const, becomes small the apex of the parabola moves far 
away and when the direction coincides with that of the ciu-ve 
a = const, the parabola degenerates into two parallel lines per- 
pendicular to the direction of the curve a = const. In this case 
the best position for the axis of x is in the direction of the curve 
a = const. Without going into any detailed investigation about 
the best position of the axis of x we can establish the general rule 
not to make the axis of x perpendicular to the direction of the 
ciu've satisfying the differential equation, that is to say, not to 
make it parallel to the axis of the parabola. But we hardly need 
pronounce this rule. In practice it would enforce its own observ- 
ance, because for that position of the axis of x not only df/dy but 
also f(x, y) are infinite and it would become impossible to plot 
the curve Y = f(x, yi). 

There is another graphical method of integrating a differential 
equation of the first order 

which in some cases may well compete with the first method. 
Like the first it is the analogue of a certain numerical method. 
The numerical method starts from given values x, y and cal- 
culates the change of y corresponding to a certain small change 
of X. Let h be the change of x and k the change of y, so that 
x+ h,y + k are the coordinates of a point on the curve satisfying 
the differential equation and passing through the point x, y. k is 
calculated in the following manner. We calculate in succession 
foiu* values kiy fe, kz, k^ by the following equations — 
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K=fix+h,y+h)h. 

We then form the arithmetical mean^ 

t. + is , h + h 

P = — 2— M-lI ! = — 2~' 

and find with a high degree of approximation as long as A is 
not too large 

^ = P + 3(9 - P)-' 
The new values 

X= x+h, Y=y + k 
are then substituted for x and y and in the same way the coordi- 
nates of a third point axe calculated and so on. 

This calculation may be performed graphically in a profitable 
manner, if the function f(x, y) is represented in a way suited to 



Fro. 81. 
' See W. Kutta, ZeiUchrift far Maihematik and Phytik, Vol. 46, p. 443. 
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the purpose. Let us suppose a number of equidistant parallels 
to the axis of ordinates: x = Xq, x =^ Xi, x = X2, x = xz, • • •. 
Along these lines f(x, y) is a function of y. Let us lay off the 
values otf(x, y) as ordinates to the abscissa y, the axis of y being 
taken as the axis of abscissas. We thus obtain a number of 
curves representing the functions /(a^o, y), f(xi, y), f{x2, y), • • •. 
Starting from a point A(xo, yo) on the first vertical x = xq (Fig. 
91) we proceed to a point Bi on the vertical x = Xz'm the following 
way. By drawing a horizontal line through A we find the 
point A' on the curve representing f{xo, y). Its ordinate is equal 
to f(xo, yo). Projecting the point A' onto the axis of x we find -4" 
and draw the line PA'\ P is a point on the negative side of 
the y-axis and PO is equal to the unit of length by which the 
lines representing f{x, y) are measured. Thus 

OA^'IPO = Six,, yo). 

Now we draw ABi perpendicular to PA'\ so that if h and h 
denote the differences of the coordinates of A and B, we have 

h/h = OA^'/PO, 

h = f(xo, yo)h. 

From Ci the point of intersection of the line ABi and the vertical 
ar = a;i we find Ci and Ci' in the same way as we found A^ and A'^ 
from A, only that Ci is taken in the ciu^e representing the 
values of f(xi, y), and draw the line AB2 perpendicular to PCi\ 
Denoting the difference of the ordinates of A and B2 by fe we have 

h PO -jy^^ + 2'^'^2)' 

or 

h = flxo + 2 ' yo + 2j*- 

From C2 the point of intersection of the line AB2 and the 
vertical ar = ari we find in the same way a point -83 on the vertical 
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X — X2 and the difference h between the ordinate of Bz and that 
of A is 



h =/(iPo + 2'yo + -2 )*• 



From Bz we pass horizontally to Bz' on the curve representing 
fi^y) and vertically down to Bz''. The line AB4, is then drawn 
perpendicular to PBz", so that the difference k^ between the 
ordinates of B^ and A is 

h = f(xo + h,y+ kz)L 

■ 

The bisection of B2-B3 and of BiB^ gives us the points Ei and E2 
and the point B is taken between Ei and E2, so that its distance 
from El is half its distance from E2. The point B is with a high 
degree of approximation a point of the curve that passes through 
A and satisfies the differential equation. 

B is then taken as a new point of departure instead of A, and 
in this manner a series of points of the curve are found. 

In order to get an idea of the accuracy attained the distance 
of the vertical lines is altered. For instance^ we may leave out 
the verticals x = Xi and x = Xz, and reach the point on the 
vertical a: = 2:4 in one step instead of two. The error of this 
point should then be about sixteen times as large as the error on 
the same vertical reached by two steps, so that the error of the 
latter should be about one-fifteenth of the distance of the two. 
If their distance is not appreciable the smaller steps are evidently 
unnecessarily small. 

The values of f(x, y) may become so large that an incon- 
veniently small unit of length must be applied to plot them. In 
this case x and y have to change parts and the differential 
equation is written in the form 

dx 1 



dy fix, y) • 
The values of l//(ir, y) are then plotted for equidistant values of 



136 GRAPHICAL BiETHODS. 

y as ordinates to the abscissa x and the constructions are changed 
accordingly. 

§ 16. Differential Equations cf the Second and Higher Orders. — 
DifiPerential equations of the second order may be written in the 
form 

dx" 



'/(^'^'S) 



Let us introduce the radius of curvature instead of the second 
diflFerential coeflBcient. Suppose we pass along a curve that 
satisfies the equation and the direction of our motion is deter- 
mined by the angle a it forms with the positive axis of x (counted 
in the usual way from the positive axis of x through ninety 
degrees to the positive axis of y and so on), s being the length of 
the curve counted from a certain point from which we start. 

We then have 

dy ^ dx 

Consequently 

(Py _ 1 da ^ 1 da 

da? "" cos* a dx cos* a ds * 

or 

da , d^y 

-^ = cos»«^. 

dafds measures the "curvature," the rate of change of direction 
as we pass along the curve, counted positive when the change 
takes place to the side of greater values of a (if the positive axis 
of X is drawn to the right and the positive axis of y upwards a 
positive value of da/ds means that the path turns to the left). 
Let us count the radius of curvature with the same sign as da/ds 
and let us denote it by p. Then we have 

- = cos' af(x, y, tg a). 

Thus the differential equation of the second order may be said 
to give the radius of curvature as a function of x, y, a, that is to 
say, as a function of place and direction. 




X 
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Let US assume that this function of three variables is repre- 
sented by a diagram, so that the length and sign of p may quickly 
be obtained for any point and any direction. 

Starting from any given point in any given direction we can 
then approximate the curve satisfying the differential equation 
by a series of circular arcs. Let A (Fig. 92) be the starting point. 
We make MaA perpendicular to the given direction and equal to 
p in length. For positive values of p, Ma must be on the positive 
side of the given direction, for negative 
values on the negative side. Ma is 
the center of curvature for the curve 
at A. With Ma as center and MaA 
as radius we draw a circular arc AB 
and draw the line BMa* On this line 
or on its production we mark the ___ 
point Mb at a distance from B equal 
to the value of p that corresponds to 
B and to the direction in which the j^^. 92. 

circular arc reaches B. With Mb as 
center and MbB as radius we draw a circular arc BC and so on. 

The true curve changes its radius of curvature continuously, 
while our approximation changes it abruptly at the points 
A, B, C, ' " . The smaller the circular arcs the less will accu- 
rately-drawn circular arcs deviate from the curve. But it must 
be kept in mind that small errors cannot be avoided, when 
passing from one arc to the next. Hence, if the arcs are taken 
very small so that their number for a given length of curve 
increases unduly, the accuracy will not be greater than with 
somewhat longer arcs. The best length cannot well be defined 
mathematically; it must be left to the experience of the draughts- 
man. 

Some advantage may be gained by letting the centers and the 
radii of the circular arcs deviate from the stated values. The 
circular arc AB (Fig. 92) is evidently drawn with too small a 
radius because the radius of the curve increases towards 5. If 
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we had taken the radius equal to MiJS it would have been too 
large. A better approximation b evidently obtained by making 
the radius of the first circular arc equal to the mean of MaA and 
MbB, and the direction with which it reaches B will also be closer 
to the right direction. 

To facilitate the plotting an instrument may be used consisting 
of a flat ruler with a hole on one end for a pencil or a capillary 
tube or any other device for tracing a line. A straight line 
with a scale is marked along the middle of the ruler and a little 
tripod of sewing needles is placed with one foot on the line and 
two feet on the paper. Thus the pencil traces a circular arc. 
When the radius is changed, the ruler is held in its position by 
pressing it against the paper until the tripod is moved to a new 
position. By this device the pencil must continue its path in 
exactly the same direction, while with the use of ordinary com- 
passes it is not easy to avoid a slight break in the curve at the 
joint of two circular arcs. 

Another method consists in a generalization of the method 
for the graphical solution of a differential equation of the first 
order. 

A differential equation of the second order 



^^i^^^y^t) 



may be written in the form of two simultaneous equations of the 

first order: 

dy ^ 

dz 

^ = 9(x, y, z). 

Let us consider the more general form, in which the differential 
coefficients of two functions y, z ot x are given as functions of 

X, y, «: 
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We may interpret x,p,za3 the coordinates of a point in space 
and the differential equation as a law establishing a certain 
direction or the opposite at every point in space where /{x, y, z) 
and g(x, y, z) are defined. A curve in space satisfies the dif- 
ferential equation, when it never deviates from the prescribed 
direction. Its projection in the xy plane represents the function 
y and its projection in the xz plane represents the function z. 

Let lis represent y and z as ordinates and x as abscissa in the 
same plane with the same system of coordinates. Any point in 



space is represented by two points with the same abscissa. The 
functions f{x, y, z) and g(x, y, a) we suppose to be given either 
by digrams or by certun methods of construction or calculation. 
For any point that we have to deal with, the values of /{ar, y, z) 
and gix, y, a) are plotted as ordiuates to the abscissa x, but for 
clearness sake not in the same system of coordinates as y and z, 
but in another system with the same axis of ordinates and an 
axis of X parallel to the first and removed far enough so that the 
drawings in the two systems do not interfere with one another. 
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Starting from a certain point P{xp, yp, Zp) in space we represent 
it by the two points Pi{xp, yp) and P^ixp^ Zp) in the first system 
and the values o{f(xp, yp, Zp) and g{xp, yp, Zp) by the two points 
Ax and At in the second system of coordinates (Fig. 93). The 
points Ai and A% determine certain directions MAi, and MA^ 
of the curves x, y and x, z, the point M (Fig. 93) being placed at a 
distance from the axis of ordinates equal to the unit of length by 
which the ordinates representing f{x, y, z) and g(x, y, z) are 
measured. Through Pi and P2 we draw parallels to MAi and 
MA2' as far as Qi and Q2 with the coordinates Xq, yq and Xq, Zq. 
With these coordinates the values f{xq,, yq, Zq) and g{xq, yq, Zq) 
are determined, which we represent by the ordinates of the 
points Bi, B2. These points again determine certain directions 
parallel to which the lines Qii2i and Q2R2 are drawn, etc. In this 
manner we find first approximations yi and Zi for the functions 
y and z and corresponding to these approximations we find 
curves representing /(a;, yi, Zi) and g{x, yi, Zi). These curves are 
now integrated graphically, the integral curve of f{x, yi, Zi) 
beginning at Pi and the integral curve of g{x, yi, Zi) at P2 and 
lead to second approximations ^2 and %: 

Jfe = yp + I f(^f Vh zi)dx, 

22 = Zp + I gi'^, yu zi)dx. 

For these second approximations the values of f(x, ^, 2^) and 
g(x, y2, Z2) are determined at a number of points along the curves 
X, y2 and x, zi suflBciently close to construct the curves representing 
f{x, 2/2, 22) and g{x, y2, Z2). By their integration a third approxi- 
mation yz, Zi is obtained 

ys = yp + I f{x, y2, Z2)dx, 



28 ^ ^p'^ \ 9(^' ya, ^)dx, 



DIFFERENTIAL AND INTEGRAL CALCULUS. 141 

and so on as long as a deviation of an approximation from the 
one before can still be .detected. As soon as there is no deviation 
for a certain distance x -— Xp the curve represents the true solu- 
tion (as far as the accuracy of the drawing goes). The curve is 
continued by taking its last point as a new starting point for a 
similar operation. 

The distance over which the integral is taken can in general 
not surpass a certain limit where the convergence of the approxi- 
mations ceases. But we are free to make it as small as we please 
and accordingly increase the number of operations to reach a 
given distance. It is evidently not economical to make it 'too 
small. On the contrary, we shall choose it as large as possible 
without unduly increasing the number of approximations. 

In the case of a differential equation 

d3? 



= K*'^'I) 



we have f{x, y, z) == z, and the curve z, x is identical with the 
curve representing the values of f(x, y, z). We shall therefore 
draw it only once. 

The proof of the convergence of the approximations is almost 
the same as in the case of the differential equation of the first 
order. 

For the n + 1"* approximation we have 

2/n+i = yp+ I f(x, yn, Zn)dx; Zn+1 = Zp + 1 9(x, yn, Zn)dx. 

For the true curve that passes through the point Xp, yp, Zp we 
find by integration 

f(x, y, z)dx; z= Zp+ I g{x, y, z)dx; 
heuce 

[f(x, yn, Zn) - fix, y, z)]dx; 

Zn+i - z = J [g(x, yn, Zn) - g(x, y, z)]dx. 
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Now let US write 

fix, yn, 2„) - fix, y, 2) = _ iyn - y) 

yn y 

, fjxy y, Zn) - fix, y, z) 

T ;; — I (2n - 2), 

and similariy 
, ^ f ^ g(a:, yn^ Zn) — gix, y, Zn) , . 

gix, y», Zn) - gix, y, 2) = iyn - y) 

yn y 

, gfe y, Zn) - gix, y, 2) 
i ; ; (2» — 2;. 

2n ^ 2 

The quotients of differences 

fJXy yn, Zn) - fix, y, Zn) 

Vn-y 

ana the three others are equal to certain values of cjf/dy, dfjdz, 
dg/dy, dg/dz for values of y, 2 between y and yn and between 2 
and 2» iy, yn, z, Zn not excluded). Let us assume that for the 
region of all the values of x, y, 2 concerned the absolute value of 
df/dy and df/dz, is not greater than Mi, and that of dg/dy and 
dg/dz not greater than M2, and that dn, €» denote the maximum 
of the absolute values of y — yn and z — Zn in the interval 
Xp to X. Then it follows that the absolute values of 

fix, yn, Zn) - fix, y, 2) and gix, yn, Zn) - gix, y, 2) 

are not greater than 

Miidn 4- €n) and 3f2(5n + 6n). 

Hence for the maximum values of y^+i — y and 2»+i — 2, which 
are denoted by d„^i and €,^^.l we obtain the limits 

5n+l ^ Miidn + €n)\X'- Xp\, €^i ^ 3/2(5n + 6n) | « — a^p | , 

and 

5n+l + €n+l ^ (^1 + M2)\X-' Xp\ iSn + €n). 

If therefore the interval x -- Xp of the integration is so far 

reduced that 

iMi+ M2)\x-Xp\^k<l, 
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5n+i + €n+i is not larger than the fraction k of (5n + ^n), but 
from the same reason 

(5n + €n) ^ A:(5,^i + €nr.i), (5„_i + €f*-i) < *(5»-2 + €11^2), ctc; 

therefore 

6h-i+€„+i ^ i"(6i + ci). 

That is to say, for a suflSciently large value of n S^+i and €„+! 
will both become as small as we please. 

As in the case of the differential equation of the first order it is 
not worth while, as a rule, to investigate the convergence for the 
purpose of finding a suflSciently close approximation by graphical 
methods. It is better at once to tackle the task of drawing the 
approximations and to repeat the operations imtil no further 
improvement is obtained. The curve will then satisfy the 
differential equation as far as the graphical methods allow it 
to be recognized. 

When the values of f{x, y, z) or g(x, y, z) become too large 
we can have recourse to the same device that we found useful 
with the differential equation of the first order. Instead of x, 
one of the other two variables y or z may be considered as inde- 
pendent, so that the equations take the form 

dx 1 dz __ g(x, y, z) 

dy " fix, y, z)' dy^ fix, y, z) * 
or 

<fa^ 1 dy _ fix, y, 2) 

dz gix, y,z)' dz gix, y, z) ' 

or we may introduce a new system of coordinates a:', y\ z' and 
consider the resulting differential equations. 

The second method for the integration of differential equations 
of the first order can also be generalized to include the second 
order. Let us again consider the more general case 

^=My,2), j^ = gix,y,z). 
Starting from a point x, y, z the changes of y and z (denoted by 
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h and J) can be calculated for a small change A of a; by the fol- 
lowing formulas analogous to those used for one differential 
equation of the first order: 

h = /(ar, y, z)h\ k = g{x, y, z)h; 

h = }{x + A, y+ is, 2 + h)h\ k = g{x+ h, y+ fe, 2 + fe)A; 



^2 4" *^3 h + h , k + k , h+ I 



f 



4 



and with a high degree of approximation, 

}c = v + Uq-p): i = p' + W-p')- 

These calculations may be performed graphically. For this 
purpose the functions f(x, y, z) and g(x, y, z) must be given in 
some handy form. We notice that in our formulas the first 
argument assumes the values x, x-\- h/2, x -\- h. In the next 
step where x+h,y'{-k,z-\-l are the coordinates of the starting 
point that play the same part that x, y, z played in the first 
step, we are free to make the change of the first argument the 
same as in the first step, so that in the formulas of the second 
step it assumes the values x+ h, x+ ^h, x+ 2h and so on for 
the following steps. All the values of the first argument can 
thus be assumed equidistant. Let us denote these equidistant 
values by 

Xo, Xi, X2, Xz, • • • . 

The values of f(x, y, z) and g{xy y, z) appear in all our formulas 
only for the constant values 

X = Xo, Xi, 0^2, ' - ' . 

For each of these constants / and g are functions of two inde- 
pendent variables and as such may be represented graphically 
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by drawings giving the curves / == const, and g = const., each 
value of X corresponding to a separate drawing. These drawings 
we must consider as the graphical form in which the differential 
equations are given. It may of course sometimes be very tire- 
some to translate the analytical form of a differential equation 
into a graphical form, but this trouble ought not to be laid to 
the account of the graphical method. 

The method now is similar to that used for the differential 
equation of the first order, y and z are plotted as ordinates in 
the same system in which z is the abscissa. Equidistant parallels 
to the axis of ordinates are drawn 

X =^ xq, X = xi, X ^ Xgy etc. 

On the first a; = aro we mark two points with ordinates y^ and zo^ 
and from the drawing that gives the values of /(a;oi y^ z) and 




Fio. 04. 

g{^^y Ify a) ^ functions of y and z we read the values f(xo, yo, afo) 
and g(xo9 y^i, zq) and draw the lines from xq, yo, and xq, zq to the 

points 

Zi, yo + ki and xz, zo + h. 

The intersections of these lines with the parallel x^ xi furnishes 

the points 
11 
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^u yo+ '2 ^^^ a^i> «o + 2 • 
With these ordinates we find from the second drawing the values 

f{xi,yo + -^,Zo + ^) and glxuyo+'^,zo + ^J, 

and by their help we can draw the lines from xo, yo and xq^ zq 
to the points 

Xi,yo + h and arj, yo + k- 

The intersections of these lines with the line x = Xi furnishes the 
points 

xi, yo+ ~2 and xi, zo + 2 > 
and with these ordinates we find the values 

/ {xu yo + ^, 2o H- 2 ) ' ^ ( ^1* yo + y , 2o + 2 } ' 

which enable us to draw the lines from xo, yo and Xo, zo to X2, 
yo + kz and X2, zo + k- 

With these two ordinates we find from the third diagram (x «= X2) 
the values 

f(x2, yo + h, Zo + k) and g(x2, yo + h, 20 + k), 

which finally enable us to draw the lines from Xoyo and o^qZo to 
X2, yo + h and a^, zo + k- 

On the vertical line a: =* a^ we thus obtain four points, Bi, B2, 
Bz, Bi, corresponding to yo + K yo + k, yo+ h, yo + h and 
four points, Bi, B2', Bz% 5/, corresponding to Zo + h, Zo + fe, 
Zo +kzo + h (Pig. 94). 

B%Bz and B\Ba are bisected by the points Ci and Ci; B%Bz' 
and Bi'B^' by the points d', C2'. FinaQy C1C2 and Ci'Ci' are 
divided into three equal parts and the points B and B' are found 
in the dividing points nearest to Ci and Ci'. 

The same construction is then repeated with B and jB' as 
starting points and furnishes two new points on the vertical 
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X = Xi and so on. To test the accuracy the construction is 
repeated with intervals of x of double the size. The difference 
in the values of y and of z found for x = Xi enables us to estimate 
the errors of the first construction — ^they are about one-fifteenth 
of the observed differences. 

Both methods are without difficulty generalized for the integra- 
tion of differential equations of any order. We can write a 
differential equation of the nth order in the form 





(ft»"" 


J\t,x 


'dt'- 


'' dt^-^r 


or in 


the form of n simultaneous equations of the first order 




dx 
di^ 


ih, 








dxi 
dt^ 


«», 








dXnS 

dt " 


Xn-U 






dx^i 


fit, X, 


Xu ^, 


• • • Xn-l). 



A more general and more symmetrical form is 

dx ^ 

•^ = /aC^, X, Xu • • • Xn^i), 



= fn(tf X, xi, • • • a?,^-i). 



di 

The functions x, X\, a^, • • • Xn^-i are then represented as ordinates 
to the abscissa i^ so that we have n different curves. When the 
function f(ff x, xi, o^, • • * Xfi^i) is given in a handy form, so that 
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its value may be quickly found for any given values of t, x, xi, 
• • • Xf^u there is no difficulty in constructing n curves whose 
ordinates represent the functions x, xt, x%, ••• a^n^i. Starting 
from given values of t, x, xi, Xi, •** x».i we have only to apply 
the same methods that have been explained for the first and the 
second order. 
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